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ABSTRACT 

Or. August 20 , 1978, the Synthetic Aperture Radar (SAR) 

cr board the satellite SEASAT, gave evidence of high energy 
internal oceaiic activity at the shelfbreak in the northern 
part of the Bay of Biscay. 

Quantitative spatial measurements of internal wave 
patterns were correlated with conventional, guasi-svncptic 
in situ da~a and yielded phase speeds of 0.55h/S for high 
frequency, nonlinear internal waves, generated in groups at 
the canyons indenting the shelfbreak and at the time of low 
tide in Brest. Their periods were found to be of 70MIN and 
their amplitudes estimated to be 15M. Ocean swell refrac- 
tion, observed on the SAR image, together with a localized 
significant increase in wave height detected by the SEASAT 
altimeter, was related to the shear of a northwestward gecs- 
trophic jet (ca. 0.5M/S) at break. 

Possible operational consequences of these features are 
summarized as an assessment of a SAP's capability ^o give 
tactical as well as scientific real-time information on the 
internal ocean. 
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RESUME 



Ls 20 Acut 1978, Is 5AR (Synthetic Aperture Radar) 
embargue sur le satellite SEASAT, a mis an evidence une 
intense activite interns ds 1* ocean au voisinage du talus 
continental dans la partie nord du Golfs de Gascogne. 

Des mesures qaantitativ as dss manifestations de surface 
des ondes internes ' ont ete corralees avec des donnees 
conventionnelles "in-situ", et quasi simultanees. Des 
vitesses de phase de 0,55 metres par seconds ont ste mises 
en evidence pour ies trains a* ondes internes non-lineaires, 
gsneres a l'heure de la basse mer de Brest, dans les canyons 
bordant le plateau continental; la periode de ces ondes sst 
de 7 0 minutes et leur amplitude a ate estimee a 15 metres. 
La refraction de la houle observaa sur 1* image SAR, en 
conjonction avec 1 1 observa tion par l’altimetre de SEASAT 
d'une amplification localisee mais significative de la 
hauteur de vague, a ete attribute au courant de cisaillem ant 
dun "jet'* geostr ophi gue da surface (da l'ordre de 0,5 metres 
par seconds) portant vers le nord-ouest le long du talus. 

Les consequences cperat ionelles possibles de ces phenc- 
menas sont reca pit ulees, pour damontrer les capacites d'un 
radar SAR a donner des informations ractiques aussi bien qua 
scientif igu es, en temps real, sur 1' ocean interne. 
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I. IN TRODCJCTIDN 



A. FORE5JARD 

Remote sensing gives satellites an essential role in the 
monitoring of ocean variability. Synthetic Aperture Radar 
(SAR) in particular, has the ability "to display images of 
surface manifestations of infernal processes which could 
typically only be inferred from lengthy and costly in-situ 
sampling from ships at sea or instrumented buoys. Even 
though the dynamical nonlinear linkage involved in "he SAR 
imaging mechanism is presently incompletely understood, the 
vivid and spectacular images ; e. g . , as those acquired by 
SEASAT, give a tremendous amount of information about v=ll- 
knovn or totally new oceanic phenomena. To study such an 
image of the 3ay of Biscay (Fig. 1. 1 and 1.2) was a chal- 
lenging and exciting topic of thesis research for the author 
who has experienced, over several years, the sea-going 
approach to military oceanography in this area. The excep- 
tional opportunity for doing such research was offered at 
the Naval Postgraduate School where direct suoporf was found 
and the introduction to other facilities was given, <=spe- 
cially to the Fleet Numerical Oceanography Center in 
Monterey and to the Jet Propulsion Laboratory in Pasadena. 

B. DATA ACQOISITION 

The optically processed SAR image was the starting point 
of the total data set collected. It was forwarded by th a 
Jet Propulsion Laboratory (JPL) of the California Institute 
of Technology in Pasadena, and a first heuristic study of 
the imag= gave the necessary evidence of its research possi- 
bilities. The swath is 1 OOKM wide and goes roughly from 



14 



Brest in Britanny to La Coruna la Galicia; internal wave 
activity is clearly visible, particularly over the ccr.ti- 
nental slope (Pig. 1.2). 

Conventional, quasi-synoptic data were then sought and 
requested, either directly througa the French Navy authori- 
ties, or by th s Naval Postgraduate School from U.S. and 
foreign, military and civilian institutions. 

A fairly danse set of XBT's (bathythermographs) was 
collected from the Fleet Numerical Oceanography Center 
(FNOC) and the French "Servi.ce Hvdr cgraphique et 
Oceanographique de la Marine" (SHOM) . The format of these 
BT’s gave rhe depth-temperature pairs an turning points; 
only five profiles extend below 50DM. Sixty nine X3T 
profiles were finally considered, corresponding to a time 
period of seven days around the SAR date and located in the 
neighbourhood of the image. The densest spacing covers the 
continental slope, bur not where the most intense activity 
is noted from the SAR image. 

Local atmospheric data were obtained from the Fr-nch 
National Weather Service. Global wearhsr chares over the 
Atlantic Ocean and swell hindcast for the Bay of Biscay were 
obtained from the FNOC archive center. 

Contemporaneous moored current meter time series were 
considered as important auxiliary information regarding the 
internal barorlinic modes of vertical movements, and the 
mean flow. The closest mooring was a subsurface array at 
1524, 3032, 3976M depth anchored at 46°57 , 2N, 009°57«1W 

during the Northeastern Atlantic Deep Sea (NEADS) mission. 
The digital data (current speed, direction, temperature, 
recorded depth) were forwarded ov the "Centre National 
d' Exploitation des Oceans" in Brest, where the " Bureau 
National des Donnees Ocsaniques" is 
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located 



Other satellite sensors carried by SEASAT acquired 
substantial information on the same area at rapid intervals. 
It was another privilege to get direct computer access to 
the Pilot Ocean Data System (PODS) of JPL, thus enabling a 
thorough survey of related data given by the SEASAT altim- 
eter and scatteromet er. 

Climatological data were blended from several sources in 
the open literature, recent papers and official unclassified 
documents of the CJ.S. and French Navies. 

C. OBJECTIVES 

The main aim in this study is to establish a definite 
correlation between the SAR surface patterns and internal 
quasi-synoptic features established by other means: conven- 

tional hydrographic observations, current observations and 
other satellite sensors. 

Avoiding the complex problem of physical linkage in the 
SAR imaging process, a pragmatic, heuristic approach -o SAP. 
analysis was made, thus tentatively establishing the power 
of SAR to monitor the oceans. The study has beer, restricted 
to the upper layer of the Bay of Biscay on each side of the 
continental slope and shelfbreak, where both the SAR image 
and other data bases were found to provide rich information. 
3eing interested only in internal oceanic features, the 
swell was only considered in regards to interaction with 
quasi-steady currents. 

In the course of the research, many different topics had 
to be considered, thus greatly enhancing the personal 
benefit that could be gathered strictly from the SAR study. 
Conversely, this fact, added to the limits in time and 
scientific expertise of the author, prevents -his study from 
b^ir.g in any way exhaustive. Furthermore, an inherent 
feature of this kind of quasi- s vno otic study is its 
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ALL the 



inability to support ob j set i vs "sea-truth". ALL ths results 
presented h-re must be thought of as describing a par-icular 
event at a particular time ar.i in a particular environment 
that can only be considered as part of the ever-growing 
information available to oceanograpaers. 

Overall, the objective of a ohasis is no increase the 
knowledge of the author himself and, in this case, in topics 
related both to the use of conventional methods and cf prom- 
ising satellite-related techniques. Consequently, this 
research was a learning sxDerier.ee and an enrich ire intro* 
duction to the new era cf operational military satellite 
oceanography. 



D. UNITS NOTATIONS CONVENTIONS 

Unless explicitly stated otherwise, the following units 
will be used, with the S.I. system as a general basis. 

Units of time 

GMT Greenwich time (Z) in Hours, Minutes, Seconds. 

(The local summer time in Prance is: GMT-020000; 
the solar time at the SAR image location is: GMT-002009) 
Day-Hour groups are expressed as DDHHMM , implicitly days of 
August 1978. 

Units of length 
Depths in mete rs (M) 

Distances in kilometers (KM) 

2^ t em peratura . 

Degrees Celsius (°C) . 

G eographic coordinates 

Geographic units are expressed in decimal degrees of lati- 
tude and longitude (lat,lon). Longitudes are counted from 0 
to 360 in most computations. 
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A uxiliary ge og ra phi c g^id 

An auxiliary Cartesian grid was adjusted over a map projec- 
tion of the area to meet the needs of the particular 3AR 
image orientation and general bathymetric gradient (Fig. 
1.1) This grid is in units of kilometers. The x-axis is 
oriented along the shelfbreak with positive values towards 
the southeast (directed seaward in the Bay). 

The y-axis is oriented along the center of the SAR 
swath and towards the northeast (directed onto the French 
continental shelf) . 

Hence x values give the relative positions on each side of 
the SAR image :(|x|< 60 for locations on the ground track), 
v values give relative positions on each side of the shelf- 
break. 

The transformation was established as follows : 

Two geographic points clearly visible on the SAR imacr= and a 



representative offshore 


bathymetric s 


cun ding were chosen t 


get the linear transform 


pa rameters. 




lat 47.30 


x=0 00 . 00 




Ion 355.67 


y= 1 32 . 42 


Penmate* h point. 


la* 43.59 


X=000 . 00 




Ion 352.74 


Y=- 38 8 . 9 2 


Los Cairos point. 


la* 46.00 


X 

II 

Oh 

• 




Ion 356.00 


1=000 . 00 


Offshore reference. 


The corresponding relati 


on ships are : 




la* =-0.006261 


X + 0.003075 


Y + 46.73 


Ion = 0.009205 


X * 0.005620 


Y + 354. 92 



X =-51.315 lat + 73.73 Ion - 23772. 
y = 34.050 lat + 57.17 lor. - 24217. 
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Cot. ferial biconic Drojection. 

Circles locate X3T* so i ndiras . 

Isobaths (in asters) are sracctnad to 10KS. 




GEOGRAPHIC GRID. 
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Figure 1. 1 




PHOTO REPRODtJCTIO?? OF SHE SIR IHAG2. 
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II. S A R IMAGING OF INTERNAL WAVES 



Starting with a survey of the recent 
SAR principles and its capabilities 
features, the various assessments made of 
ar.isms of interaction between an internal 
short wind-driven capillary/ultragravitv 
can detect are reviewed. 



literature on the 
to detect oeanic 
the dominant mech- 
wave field and the 
waves that a SA? 



A. SYNTHETIC ANTENNAS 



Synthesizing a radar antenna consists of using a moving 
antenna of limited size and combining all information gath- 
ered on a given point of the earth during a given time 
interval (Fig. 2.2) . The distance of travel on the orbit 
gives the length of a virtual antenna of much finer resolu- 
tion. The resulting resolution is independent of range and 
speed, and it is theoretically half the real antenna's aper- 
ture. In ~he range (cross -velocity ) direction, the length 
of coherent pulses determine (in tie usual wav) the corre- 
sponding accuracy. 

In the case of SEASAT, the data were transmitted in real 
time by the satellite to one of the special oro cessing earth 
stations. The final record is either a computer-compatible 
tape format, or a film conversion. The photographic output 
is displayed in four parallel subswaths to accoun* for a 
progressive decrease in the radar response intensity with 
range (Fia. 1.2) . Some particular images over limi~ed areas 
have been digitally processed with a valuable enhancement of 
resc lut ion. 
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The extraction of a t w o-dimensio nai surface linage spec- 
trum is done either by an optical Fast Fourier Transform 
[Gonzales 1979 and Shuchman 1978], or by numerical methods 
[algorithm of Hasselman 1990]. Considerations of the rela- 
tive qualities of these processing techniques can be found 
in Mitchel [ 1992 ] and Tomiyasu [ 1 978 ]. 

Characteristics of the SEASAT SAR are given here {Table 
I an d Fig. 2. 1 ) . 



TABLE I 

SEASAT SAR CHARACTERISTICS 



| satellite altitude 


| 3 03 


KM 




j satellite velocity 


| 7.2 


KM/S 




j near-polar orbit 


| 1 08 


degre 


33 | 


j view angle from nadir 


I 20 


degre 


3 3 


I swath width 


I 1 00 


KM 




| wavelength 


.235 


M 




| frequency L-band 


| 1275 


GHZ 




1 real antennae 


| 1 1 x2 


M 




svnthetic antenna 


1 15 


KM 




I antenna polarization 


| 9 9 % 


1 - r> oa 


r he r i z • j 


j centra st ' ratio 


1 l 


d 3 ~ 




1 radar cain 


1 35 


d3 




j processing resolution 


1 25 


M 




I bulse repetition rate 


I 5 00 


mi cro 


S9eo r.ds | 


1 bu l=s length 


I 33.4 


micro 


s 9 z o n d s j 



B. OCEAN IMAGING MECHANISMS 

The application of SAR systems to oceano err a phy requires 
extension of the basic SAR earth imaging theory to account 
for scattering from a time-varying surface that satisfies 
nonlinear hydrodynamic constraints. On the other hand, the 
general scattering models for low resolution radars do not 
apply to the spaceborne high-r esoluticn SAR's. Th= 
complexity of the problem results from the widely nonlinear 
electromagnetic and hydrodynamic interactions. But this 
nonlinearity is at the same time one main reason for the 



remarkable capabilities of SAR far detection of mesoscale 
oceanic phenomena. 

1 . Dopo le r Effects from the Roving Ocean Surface 
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optimal focus of SAR imagery over ocean surfaces, the b a st 
speed correction to use is the long nave phase velocity, and 
not as would seem intuitively, the long waves orbital 
velocity on which the capillary shore waves ride. In the 
case of internal waves, this phase velocity is given by 



Cp =gh 



ft - p. 



y 



? J e* 

where "h” is the depth of the thermocline and "P,", "Pa.” are 

the densities each side of the thermocline. This displace- 
ment effect will periodically modulate the image intensity 



long wave pattern. In the case of 
(cross-track component toward the satellite) the 
image will be artificially translated 
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tion will increase accordingly and there will be stretching 
of the image and a lower 3r aaa-sca: t er than expected. 
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2. Scattsr in.g of Radar Signal from Began tfayas 

Bas-d on a r=vi=w [Valar.zusLa 1978] of site - romag- 
r.atic wavs scattering from ocean waves, the standard 
approach considers the ocean sarface as statistically 
described, but a deterministic formulation might be more 
appropriate when the resolution area is small enough to 
correspond ~o a few wavelengths only. The former d<=-scrio- 
tion is nevertheless easier and most of the rims more accu- 
rate. The statistical approaci may lead either to a 
3raag-rescnar.oe-like mechanism or a "facet" mechanism 
(tangent- plane theory) [Brown 1976 ]. 

2* Bragg Nodal ver s us Tar, cent Plane Mode], 



The 3ragg mechanism is limited to rough surfaces 
satisfying the Rayleigh criterion : 



2. Ic W (rm s') CoS 8 4: 11 y 

where "H" is the typical surface height fluctuation and 
is the incident angle cf the radar wave on the water 
surface. The total backscatter is then given by : 



cr 



(9,1c). . = Uirk^co^S I 3 .. (9)| 



y 



where is the two-dimensional Four! 

th a sea surf as 9 with an oca an wave n amber 11 2 k 
represents ^ he radar reflectivity of the a 
face. The alternative approach is the ta 
facet model applicable to gently mdulating 
and small (Q<< 30°) incidence angles (for 
230) . 

We have in this case 



-r spectrum or 
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where R ( 0) is the Fresnel coefficient for normal incifsr.es 
ar.d "<s 2 >" is the total slops variance for ocean waves of 
long period. In both cases J *he intensity variations of the 
SAR signal will be correlated to ins energy density of 

the ocean waves. (They will be ar the same time modified by 
the mean incidence angle, resulting in a decreasing average 
intensity along a cross track line.) Therefore, the nonli- 
near periodic energy tranfer from the capillary w a v ® s to the 
mesoscale phenomena (currents, swell and internal waves) 
will dramatically affect the S AR imagery. 

C. COUPLING WITH INTERNAL OCEAN DTSANICS 

Small coherent disturbances in tne current field can 
cause important modulations in the ultragravity wavefield 
[Mollc-Christensen 1982] : e.g., for wavelengths of 30 ems, 

a shear current of 0.003 M/s ec can induce a modulation at 
the surface to which a S AR will be sensitive. 

1 . The Sur f ace Str a ini no and tne lunching Effects 

The possible mechanism involved appear to be either 
an hydrodynamic effect (the surface straining) , or a Doppler 
effect (the bunching) . 

The straining effect applies to all internal ocean 
phenomena (internal waves, shear currents, eddies). 
Variations of " Q -0 " along the direction of wave propagation 
are induced by hydrodynamic interactions between short reso- 
nant waves (Xt£.3M) and large ocean periodic motions 
(X>100M). This mechanism can be thought of as an alter- 
nating pat -or n of divergence-convergence of the free 

surface, riding on a long-wave modulation ( sw -11 or internal 
wave) . The resonant waves are saturated locally at the 
region of maximum convergence of velocity, ahead of the wave 
crest (Fig. 2. 4) . 
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In the short wav? (Phillipsi spectrum, the modula- 



tion transfer function can be 'written. [ A Ipers , Ross and 
Buff enach 1 981 ]: 




where "H” is a time constant related to the resonant waves. 



ficial displacement of the image of periodically moving 
surface [Alpers and Bufenach 1979, Alpers ei al. 1981]. 

Simply put, the Doppler shift due to surface velocity is 
attributed to location (azimuth) and the sian.al is rnislo- 

cated in the image. It seems to be the best model for 

azimuth traveling swell and eventually for internal waves. 
However, it is a geometric effect that does non account for 
energy transfer as well as 'he surface straining dynamic 
me chan is m. 

2. A oo 1 ic a t ion to Internal Waves 



may therefore be written more generally as ] Vesecky and 
Stewart 198 2 ]: 



the flight (azimuth) direction and 
the direction of propagation of the internal wave, and "Jl" 
is the angular frequency of the wavs component. 




The velocity bunching effect is she result of arti- 





wit h 




where " B " is the total modulational transfer function and 
can be considered as the result of • wo different mechanisms 
(neglecting the tilting effect appLicaole only to swell) : 



B = Bfnydrodyn.) + 3 (vel. bunching) 
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wave 



With a ^hecrv relating the surface wave field to the 
internal wave amplitude and surface winds, information on 
the internal wave amplitude, and simultaneously on the wind 
field, can be obtained from the imagery. However, such a 
theory is not available yet [Beal 1991]. 

It is generally believed that internal waves are 
seen by a radar as a result of two mechanisms : 

(1) The surface currents associated with the internal wave 
field sweeD together surface oils and materials in the 
surface water convergence zones resulting in severe damping 
of short capillary/gravity waves to form slicks, and 

(2) The surface current field induced by the internal waves 
result in a periodic modulation of the height of the 
capillary/gravity waves. 

Actually the two possible mechanisms implied give 
contradictory results: damping oy oil slicks will occur 

over the convergence lines, whereas the current hydrodynamic 
damping will be localized on the divergence areas. It is 
the author's assumption that whenever the surface patterns 
visible on the S AR image are linear, localized and filamen- 
tary (Cf patterns described in 7-3-1) , the oil slick expla- 
nation is appropriate, whereas alternation of high and low 
brightness can be mere easily explained by i y drodyr. ami cal 
damping due to internal waves (Cf group pat - erns described 
in 7-3-3 and 7-B-4) . 

According to the surface straining mechanism, an 
internal wave traveling at phase speed " c" produces a 
surface straining and a velocity pattern "u(x-ct)" moving 
with the waves. Short surface gravity waves travelincr in 
the same aeneral direction as the internal waves car. 
strongly interact, particularly when the group velocity of 
the surface wave is close to the ohase velocity of th- 
irternal wave. 



27 



o _ 



Gargett and Hughes [1981] give another de scrip- ion 
f he same principle. The horizontal surface currents 
associated with internal waves are Large, being of the same 
order as the wave phase velocities themselves. They are 
oriented approximately either parallel (in the troughs) 
antiparallel (on the crests), to the direction of propaga- 
tion. A wind f ield ' super! mpcsed on the current will excite 
a capillary/uLtragravity wave spectrum that is steeper where 
the wind opposes the internally generated current than where 
it is parallel to it (Fig. 2.U). Eventually, a modula t ior.al 
interaction couples the surface-wave spectrum to that of the 
internal waves and possibly creates a growth of the latter 
[Dysthe 1981]. 

Indirect detection of internal waves may be possible 
by analyzing the variations of the swell energy spectrum 
directly extracted by FFT of the SAR data. Refraction of the 
swell is then the mechanism by whioh th a internal waves are 
manifested on the sea surface. 

D. S 5HSITI VIT Y OF A SAR TO INTERNAL MOTIONS 
1 * 5 nv i ro n m en t a 1 Limit ati on s 

For possible radar imaging of ocean internal wave 
motion, the first essential condition is the existence of 
surface capillary or ultra -gravity waves. This excludes 
total calm wind conditions, but not rough and even breaking 
seas. The relative directions of the short surface waves 
and th- internal wave field affects the hydrcdynamical 
interaction process. Furthermore, the displacement of the 
satellite may induce Doppler imaging effects dependent on 
the angle which the radar waves make with the surface waves. 
Hence , for a given infernal wave field and known satellite 
orbit parameters, the optimal wind velocity conditions ■could 
be stated. It sterns, but if would need exoer imenf a 1 confir- 
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mation, that a wind blowing parallel to the propacaticr. of 
the internal waves is optimal for strong interaction. The 
radar modulation is largest for internal waves traveling in 
azimuthal direction, but on the other hand, it is assumed 
that a satellite track parallel to the same internal wave 
field may induce enhancing Doppler- ef facts . 

That these geophysical factors are sometimes contra- 
dictory does not seem, however, to inhibit the capability of 
a SAS to disDlay surface manifestations (observed in many 
different cases) of internal motions. In particular the 
analysis conducted in chapter V shows internal waves propa- 
gating in almost all relative directions, thus proving that 
at least seme interaction and radar modulation exists for 
every ccnfiauration. 

2. Propitious Internal Wave Prooerties 



The group velocity of wind-driver, short waves is 



given by 




9 

lie 



i<T 






where "T" is the surface tension. The optimal internal wav<= 
phase speed should then be of the same order such that 



q H 



A? j 

p v y 



implying approximately for ultra-gravity waves 

AL ~ _i_ , 

f 2.kU / 

which, in the case of the 30CH wavelength waves that S2ASAT 
SAP can detect gives 

e ~ voh 



For a two-layer ocean with a thermocline at 




s. 



5 OH we obtain 
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The surface horizontal particle vaLocities must, or. the 
other hand, ba larga enough. But we know they ara directly 
related to tha amplitude of the internal wave itself. 

We therefore think that detection of internal waves 
by their surface manifestation should ba favored by: 
a deep mixed layer (H) , 
a sharp thermocline, and 
a density jump of tha order of 1 / 43 H 
(H in meters) 

3. F avora bl a Cases : an Assessment 

Infernal tidally generated waves on the continental 
shelf have frequency properties and amplitudes that cause 
them to fall into the stated favorable conditions. Those 
internal waves are generated most commonly in the highly 
stratified seasonal thermocline by interaction between a 
barotropic internal tide and bathymetric irregularities. 
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example, the 
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offshore Fran 
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ward. The day (20 August 1978) of passage of SEASAT ever 
this area corresponds to a spring tide, and the bathythermo- 
graphs show the strong stratification common to the summer 
season. A search of meteorological data has confirmed 
average summer conditions to have prevailed, with probable 
winds from the northeast. Previoas conventional studies 
have shown that internal waves are frequently observed and 
propagate perpendicular to the 'itf/SE orientation of the 
continental slope. 
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COMBINED UTILIZATION OF CONVENTIONAL SENSORS 



Adequate knowledge of wind conditions is necessary for 
the reliability of any s AR analysis. Ship observations, as 
well as data from other remote sensors (e.g., a scatterom- 
eter) , were used in this study. (Both the scatt ercmete r and 
SAR are side-looking sensors. Hence they offer generally 
overlapping coverage.) The required knowledge of the stra- 
tification and of the spectrum of periodic movements can be 
derived from bathythermographs, and from moored instrumenta- 
tion measuring current and temperarare time-series. 

1 • B athyt he rm og rap hs 

Or.® bat hythermcgraphic or bathycelerimet ric measure- 
ment allows the computation of the eigenvalues solution to 
the internal wave generation problem 'Baines 1982]. It is 
then possible to find the wavelength of internal waves from 
the spacing of wave-packets generated by tidal flow over 
topo graphy. 

If the repetition, rate of the forcing " r " is known 
(tidal internal waves for instance) one can derive from the 
spacing between packets "a" , the pnase speed 

= cl c 

A bathythermograph in situ makes possible the computation of 
the main density gradient on the thsrmocline from 

Comparing observed and theoretical values of the phase 
speed, it is possible to infer the tnsrmocline depth by 

h - ilri . 

9 Ap y 
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and eventually, by combining this information with sea- 
surfacs temperatures, the heat storage in the ocean 
[Mollc-Christensen and dascarer.has 1 979 ]. 

2 . C ur r en t and Tempera ture Observations 

Time series of current and/or temperature measure- 
ments at several depths make possible the vertical modal 
analysis of the vertical and horizontal movements at a 
point . 

High resolution observations of the sea surface 
temperature by airborne or spaceborne infrared radiometers 
might detect horizontal gradients due to variations of the 
surface thermal layer depth directly related to the 
stretching effect of an internal wave field. 
Quasi-simuli aneous temperature transects, in the sam a way, 
can resolve the space characteristics of an internal wave 
field that can only be inferred from a fixed observation 
station together with internal wave theory. 
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^Spacecraft velocity 
x vector v~~ 8 km/s 



Synthetic aperture^ 
length r 

Z. ~ 16 km JJ?- 




Nadir 



D, real aperture, ~ 12 m 



6 = real beamwidth 

o' = effective angular 
resolution 

h , slant range, ~ 10 s m 



r a , azimuth resolution, 
~ 6 m 

= one-half real aperture 



From Spacaborr.s 3A3 far Oceanography [1981] 



Figura 2.1 S2ASAT 5A2 GEOMETRY. 
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III. GENERAL GEOP HYSI CAL FEATURES OF THE BAY OF BISCAY 

According to Sverdrup et al. '1942] "the terminal 
branches cf the Gulf Stream system are not all well known; 
irregular branches which turn to the south terminate in 
great whirls off the European coast 11 . These ir r ecu lar it ies 
could correspond to semi-permanent circulation patterns or 
propagating eddies. Recent studies have be=n made of th a 
mean circulation and of mesoscale features in the Bay cf 
3iscay, particularly in the vicinity of the Armorican conti- 
nental slope, from which some general conclusions can be 
dr awn. 

Several general geophysical features relevant to the 
subarea covered by the SAP image will be considered, giving 
particular emphasis to the summer climatological conditions 
under which the SAR image was taken. 

A. ATMOSPHERIC CIRCULATION 

Atmospheric conditions over the 3ay of Biscay are deter- 
mined by the Location of the North Atlantic Subtropical High 
(pressure center) centered in February around 25°N and in 
August around 35°N. The High extends over 20 to 25 degrees 
cf latitude in summer (vs. only 15 in winter) and is limited 
to the north by the Icelandic Low. The general atmosph-ric 
circulation is consequently a mean eastward flow, weak in 
summer, which is modified by the presence of the European 
continent and perturbed by the passage of fronts originated 
ever the Ncrt vest Atlantic in the vicinity of the Polar 
Front (between 45°and 50°N in summer). Table II aives the 
s- at ist ical averages for August at thre- stations, r e sd ac- 
tive I y north, east, and southeast of tne Bay of 3iscay. 
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BATHYMETRY 



The Eay cf Biscay (Gclfe de Gascogne) is a fairly large 
bay. It is 4500 (250NM) wide at its western opening on the 
West European 3asin and extends 5530 (330 NM) to the south- 
east (Fig. 3.1). Its physiographic characteristics hav= 
been presented by Laughton [ 1975 ]. 

The abyssal plain has an average depth of 4500M, 
increasing southward to 439 OM. It is limited to the south 
by the narrow (400) Cantabrian continental shelf, and to 
the north by a wide (2900) and shallow ( 1 3 0 M) shelf (Fig. 
3.2) . This latter shelf extends far beyond th e Bay of 
Biscay to encompass the British IsLes. The designation cf 
Celtic Sea is generally given to the northern part of the 
shelf; the Armorican Continental Shelf is to the south of 
48°N, the latitude of Brest. 

The "La Chapelle Bank", at tie above approximate limit 
(4 8 0 N , 007OW) is 100 to 1 5 0 M dee? and characterized by large 
scale and smooth ridges (amplitude ca. 40M) oriented NNE-SSW 
and 200 apart. Its southvesrard Limit (47041*8, 007oi3'W), 
sometimes referred to as the bank per se, exhibits sand 
waves [Cartwright 1958], correlated with scattering layer 
undulations [Stride 1960] as a consequence of internal tidal 
wave activity [ Carruthers 1 963 and Cartwright 196 3 ]. These 
sand-waves have an average height of 3M and a wavelength of 
859M; they are oriented toward 113°T (parallel to the edge 
of the shelf) and extend over at least 200 from the slope. 
Southeast of La Chapelle Bank, the continental seafloor has 
a variable geological nature (rocks, mud, gravel and sand) , 
homogeneous in depth with a v=ry gentle slope (3.05%) . 

The shelfbreak of the Armorican Shelf is oriented 130°T 
and occurs between 180 and 200M (Fig. 3.3) The width of the 
continental slope decreases from 350 at La Chapelle Bank to 
200 at 470, 3 05°30 , W, where the steepness is maximum (29% 
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between the 230 and 1000M isobathe') . This steep slop= is 
indented by many canyons and ridges normal to its axis. The 
designations given by Laughton [ 1 973 ] will be used here. Or. 
the swath covered by the SAR image, rive major canyons are 
visible: the Per.marc’h Canyon, me Ouilvir.ee Canyon, the 

Odet Canyon, the Blavet Canyon and the Belle-Ile Canyon 
(Fig. 3.3). 

The continental rise, south of the slope, is 200KM wide. 
It is broken by the Gascony Knoll at 45°20'N, 005°20'W, with 
a minimum depth of 3950M, at the eastern edge of the SAR 
image. The Cantabrian Knoll centered at 003°00'W 

is steeper and higher (3000 M minimum depth, 20 by 400 hori- 
zontally) , but it is 500 off the western edge of the SAR 
image . 

C. OCEAN CLIMATOLOGY 

The surface circulation in the Northeastern Atlantic is 
not as well defined as that in me western (Gulf Stream) 
regime of the North Atlantic anticyclonic circulation. After 
crossing the mid-Atlantic ridge, me Gulf Stream loses its 
coherency and splits into several o ranches, which are weak 
and greatly variable on all space and time scales and which 
eventually decay into fields of poorly defined eddies (Fig. 
3.4) . The easternmost of these branches flows south along 
the European coast, probably creating a weak cyclonic circu- 
lation in the 3ay of Biscay, as it weakens along the north 
coast of Spain [Millet 1930]. A transient westward surface 
current ir. the northern part of the Bay follows long periods 
of eastward winds; if is referred to as the Ren nel Current 
and described in the French " Instruct ions Naufiques" [SHOM 
1975] as being probably narrow, but strong (0.5 to 0.7M/S). 
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The nsar-oo ttora circulation is highly variable, t to. 
However# ■‘•.here is a quasi permanent cyclonic bottom circula- 
tion on the abyssal plain centers! at 44°30'N, 004O30*w, 

which is driven by the eastward flow along the Iberian 
Peninsula. On the Armor ic an Shelf the mean bottom current 
drifts to the northwest from March to July, and then towards 
the southeast from August to November. The average current 
speed measured at 1M above the sea floor is 0.04M/S. [Camus 
1971 ]. 

The intermediate waters in the abyssal plain of the Bay 
are under the influence of the Mediterranean outflow, which 
spreads north out of the Straits of Gibraltar. This 
intruding water mass is centered at a depth of about 1100M 
and its T-S char act eristics , originally T=11.9, s = 36.50°/oo, 
decrease by mixing to T= 10 , S=35.33°/oo in the greater part 

of the Bay of Biscay. Assuming a level of no-motion at 
2000M, results from two separate hydrographic data sets 
(July 1963 and November 1967) exhioit the following charac- 
teristics of the Mediterranean Water circulation ( Le Flcch 
1969]: a permanent cyclonic circulation occurs in the 

southwest portion of the 3a y of Biscay, approximately 
centered at 44O30'N, 007O33'W on which is superimposed vari- 
able meanders and probable eddies; the speed of this circu- 
lation is of the order of 0.0 5 M/S. 

D. WATER MASSES, FRONTAL STRUCTURES, AND EDDIES 

Hydrographicallv , the deep ocean waters of the 3av of 
Biscay have three main components; the relatively warm 
surface wafers are from the North Atlantic anticyclonic 
gyre; the deep and bottom wafers are derived from Antarctic 
Bottom and Intermediate Water; t.ne midcolumn wafer mass 
charcter i st ics vary with the degree of mixing of 

Mediterranean Waters with either one of the first two 
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'components. On the Armor! can Shelf a mass of colder, less 
saline water occars as a bottom layer adjacent to the 
French coast (Fig. 3.5). 

We will consider as a front, the locus of sharp hori- 
zontal gradients in temperature or salinity at an arbitrary 
depth. Two main surface fronts are driven by tidal 
currents, ore offshore Britanny due to a horizontal decrease 
in mixing, and the other at the continental slope which will 
be discussed in more detail. Two pelagic fronts have beer, 
identified : (1)the the rmalhaliae front delimiting the 

extension of Mediterranean Waters, which is difficult to 
monitor but which may have important acoustical effects, and 
(2) the thermal front bounding the coastal bottom cold water 
mass of the Acmorican Shelf, which is enhanced in summer by 
surface heating, but which has no surface temperature mani- 
festation in that season. 

Eddies occur in the Bay of Biscay in summer conditions, 
when the thermal structure of the upper water column is 
strongly developed. Based on infrared satellite images of 
eddies bordering the Ushant frontal zone, these eddies are 
most apparent along geostropnically balanced fronts where 
cross frontal mixing takes the form of cyclonic disturbances 
[ Pin are e 1978]. They have a length scale of 20 tc u o KM, 
with a time scale of a few days and translational velocity 
less than 2CM/S. In the case of the Armorican shelfbreak, 
irregular eddies form on each side of the 230M isobath: in 
this process "the cooler water sinks below the warmer water 
ar.d is stretched vertically" [Pingree 1978 ]. 

Over the abyssal plain, less apparent eddies have been 
detected by drifting-buoy experiments and by satellite 
imagery (Diskson 1981). These eddies, both ant icyclor. ic and 
cyclonic, seem to partly decouple salinity and temperature 
features from density surfaces. From an A'/HRS imag^ taken 
on 14 May 1979 , and in the vicinity of the SAP. image, it 
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appears that cyclonic eddies were aligned on the continental 
rise east of 006°w, and that a pain of cyclonic and ar.ticy- 
clonic eddies was located, respectively, northwest and 
southeast of 45°30'N, 007O30'W. The length scale cf these 

deep-water eddies is 50 to 100KM [Dickson 1981]. 

E. TIDES 

"The tides of the Atlantic Ocean can be considered as 
composed of two parts, a co-cscillat ing tide which is main- 
tained by the tide cf the Antarctic Ocean, and a free fide 
which is maintained by the direct effect of the tide 
producing forces. As a result, the semidiurnal fide can be 
considered as composed of two standing oscillations having a 
phase difference of one quarter of three lunar hours." 
[Sverdrup 1942] The amphiiromic point is "estimated reli- 
ably" [Cartwright 1980] at 5 0 0 M , 039OW, so that the fide 

will propagate from south to north in the Bay of 3iscay with 
a 90 degrees Greenwich phase lag in the south versus 110 in 
the north (Fig. 3.6). The tidal amplitudes increase from 
1.0H offshore, to 1 . 2X over the continental rise and 1.4M at 
50Kd off the French ccast. By way of comparison, the diurnal 
tidal amplitudes over the area are of the order of 0.07 m 
(Table III) . The tidal energy dissipation computed by 
Cartwright [1980] largly exceeds direct estimates based on 
bottom friction, leading iiu to conclude that "a substantial 
loss must occur along the Biscay shelf south of Ouessanf 
(Ushant Island) ." And indeed, ?i.igree [ 1982] infers from 
both observations and nonlinear tidal models -hat "about 503 
cf the tidal energy flux onto the shelf occurs in the 
regions where the tidal streams reach about 40C*VS (ca. 
005030' * to 309 030’ W)". Furthermore, he notes that »»the 
shelfbreak bathymetric contours tend to incr-as- this tidal 
strength wherever they are concave (canyons)". 
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TABLE III 
TIDAL CONSTITUENTS 
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Internal tidal modes have been studied from current 
meter observations at 47ou7N 010°13*S in April 1978, [James 
1982]. Results support u he notion that internal tides are 
dominated by the lowest, vertical modes and chat the M2 
component is significantly coherent with the equilibrium 
tid=. Barotropic tidal current elLipses were found to be 
narrow and directed perpendicular to the shelf, thus 
"forcing the isopvcnal surfaces over topographic features 
(the accepted explanation for the generation of internal 
tides)" [James 1992]. They show a regular turning (counter- 
clockwise) of their orientation axis from 570 to 3670M depth 
and could thus produce baroclinic activity. The tidal 
ellipses are themselves of clockwise rotation at depths less 
than 1270M. 

F. INTERNAL SJA72 MOTIONS 

Internal waves are a common phenomenon in the stratified 
ocean and exist in the frequency band from the local iner- 
tial frequency (15. 4H period at 48^ Mi to the maximum 3 runt - 
Vaisala frequency. At tidal frequencies (e.g. , such as, 
semidiurnal), '•base periodic motions should be called 
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internal (barcclinic) tiles as the/ result from forcing by 
tidal forces. At frequencies greater than th=> 3runt-7aisala 
frequency, perturbations cannot propagate as waves and are 
turbulent. 

In the center of the 3ay of Biscay, temperature time 
series f Hecht and Hughes, 1971] have shown that semidiurnal 
internal tides are greatly enhanced by spring {barotrcpic) 
tides and can generate high frequency oscillations with 
marked peaks at periods of 1=7 1 and 24HIN. These deep-water 
internal tides could be generated by the incidence of the 
barotrcpic fide on the continental shelf, to -her. propagate 
offshore normal to the coast line * Rattray I 960 ], Their 
amplitudes would then be given by: 



A 








/ 



where 



, Z [1 

l r\ - 1 

and "A" is the amplitude of the i.nt 
tidal range, ’’h" is the depth to the 
deep-water depth, "d" is the dept 
constant) , "W" is the width of the o 
the frequency of the internal wave 
parameter, "g" is the acceleration 
" h f_ s relative density dif 

layers. H=cht's [1971] observed 
estimated by the above formula ( 



ernal wav=, "2H” is the 
thermccline, "D" is the 
h of the shelf (assumed 
ontinental shelf, is 

, " f " is the Coriolis 

due to gravity, and 
ferer.ce between th® two 
amplitudes exceed those 
231 obseved vs. 15. 3h 



computed at spring fids). On the continental shelf (46°!T 
003°50W), nearbottom current records occasionally exhibit 
bursts of periodic oscillations (ca. 30MI1I) followed durina 
a few hours by apparent white noise fluctuations [ Le Finch 
1959 ]. 
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IV. GEOPHYSICAL SITUATION ON AUGUST 20 1978 



A. HETEOBOLOS Y AND WIND SI HESS 

On the 20th of' August 1978, a high pressure ridge 
extended from the Azores to Western Europe. The average 
surface pressure over th=> Bay of Biscay was 1023MB (Fig. 
4. 1) . A strong storm generated on the 15th of August 
offshore Labrador was dying over Iceland (Fig. 4.3) . A weak 
cold front was 300NM west of the Bay of Biscay on the 19th 
at 1200 GMT and probably underwent frontolysis in the area 
of interest. 

From ship reports, the evolution of the in- situ atmos- 
pheric situation (ca. 470N,005°W) can be analysed (Table 
IV). Between 200300Z and 200600Z the cold front was present 
near the location of observations at 45°40N, 004O30W. 



TABLE IV 
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From synoptic charts and these ship repor*s a comparison 
can be made to select the SEASAT sea tteromet er winds the 
time of the SAR image (fig. 4.2) . At 202142Z, a local high 
pressure center was in the vicinity of 48°N, 010° w, *hus 

establishing an anticyclonic surfaoe flow decreasing from 
ca . 10H/S south of Penmarc'h (north of the SAR image) to 

ca. 5 M/S at 45°N 006°W (south of tie SAR image) . The best 
estimate for the direction is a flow towards 200°T (winds 
north- northeast) . 

B. SWELL 

A heavy (mean significant wave height of 2.5M), and long 
(mean wavelength of 230M) swell had originated 20000 away 
(Fig. 4.4) and was currently prevailing over -he sea state. 
It was propagating roughly 130°T (from shipboard meteorolo- 
gical reports). The generation areas were estimated by 
atmospheric charts (Fig. 4.3) to be in the southwest quad- 
rant of the surface pressure low. Distances and time of 
propagation with respect to the SAR image location and time 
yielded estimates of the wavelengths that could be expected 
in the Bay of 3iscay f August 20, as a function of mean- 

direction. The deep water surface gravity waves equations 
used are the following: 



elapsed since the original wind generation. 

No swell spectra were available to compare with the 
patterns visible on the SAR image, but a correlation with 



the period is expressed as 
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thus yielding 



X = J>1 
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with "D" , distance from the gen 




and "fc '* 




altimeter pass of 21072 1Z (ascending nods passing by 46°N, 
005°W) was a sudden increase in significant wave heigh* from 
an average of 2M to more than 4.1 over a short distance 
(15 0) at position 47oi5'N, 006°00’tf (Fig. 4.5), i.e. near 
the shelfbreak. This astonishingly abrupt variation could 
hardly be related to -any atmospheric mesoscale process. It 
will be shown (Cf VI-E) that it may be the consequence of a 
strong interaction with current shear occurring at the 
shelf break. According to sailers experienced with this area 
[Chiari 1982], crossing the 2901 isobath is often accompa- 
nied by a characteristic change in sea-state. An historical 
reference is given by La Dar.cis [1943], who attributes the 
denomination of La Chapelle 3ank to an observation mad c by 
the "Sieur de La Chapelle'’ in 1695 of surface turbulence 
wrongly associated with extreme shoaling. This apparent 
shoal was resighted in 1754 and a rock was marked on the 
charts for the following decades. A well known example of 
such a current-swell interaction is the occurrence of 
"aiant" wav^s [Huthnar.ee 1 980 ] near the shelf edge off the 
East coast of South Africa. These are thought to arise at a 
caustic associated with the refraction of surface waves by 
the Agulhas Current. In :he same way, an alongshore current 
at the shelfbreak in the Celtic shelf could explain the 
increase in wave height observed there. 

C. MEAN AND TIDAL CURRENTS 

Data available for direct estimation of the mean current 
and its variations over the area, are provided by the deeo 
water subsurface mooring at 46°57 , N, 009°57* iw with current 

meters at 1524, 3032 and 39761 (3530 off the 5AR ground 

track and 1800 offshore of the shelfbreak) (Cf 1-E). Lew 
frequency currents at the three levels do not exhibit any 
obvious correlation and vary over all dir-cticns with a 
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time-scale of three months. The mean currents in the last 
10 days before the 20th of August are respectively: 
at 1524M toward OIQOT at 0.05M/S 
at 3032M toward 270®r at 0.094M/S 

(veering 33D°T in last 3 days) 
at 3976M toward 3l5or at 0.001M/S 
On a shorter time scale, progressive vector diagrams show a 
semidiurnal periodicity at each of the three levels (Fig. 
u.5) ; the m^an speeds are : 
at 15 2 ,4M 0 .0 3 .1/S 
at 30 3 2M 0.05 M/S 
12 39 76M 0.00 1 5 M/S 

At 20 21 40 , time of the SSASAT passage, the instantaneous 
current vectors were : 

at 1_5 24M toward 315°T a 4 : 0.08M/S 
11 3032M toward 18 50T at 0.06M/S 
a 4 3976M toward 270°r at 0.0015M/S 
The passage of SSASAT on 20 August followed a spring tide 
ranging 110 or. the French coefficient scale [SHOM 1968 ] 
(i.e., amplitudes 1535 greater than the average spring tides 
of coefficient 95). At the time of the imag-, the surface 
tide at the shelfbreak was at high tide plus four and a half 
hours (Table 7). The surface currents were estimated 
towards 230 or at 0.85M/S (Table 4.5|, [SHOM 1968]. 

D. SPATIAL THERMAL STROCTtJRE 

The thermal structure over the area was established from 
a set of 69 X3T f S (Cf 1-3) . The rime span of this tempera- 
ture sampling was from 0900 17 August to 0500 24 Auaust. 

The spatial Limits are between 43° and 48°North, and 10° 
West and th= coastline. The distribution of this sampling 
is conveniently centered around the radar image, in partic- 
ular in th- vicinity of the continental slope and with a 
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TABLE V 
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Days of 
Auciust 


In BREST | 
GMT Hour) 


At 

bre 


Shelf 

ak 


Coef f 


I Amo lit ud 
| in' Br est 


e f AmDl itu de 
I a 4 : ' Break 


1 5 


30 


33 


00 


13 


63 


I i . 9 3 a 


I 1.25M 




13 


39 | 


12 


49 


5 5 


| 2. 12 


1 1.38 


1 6 


3 1 


42 1 


01 


22 


73 


| 2.34 


1 1.52 




14 


12 | 


13 


52 


33 


| 2. 57 


I 1.67 


17 


32 


40 | 


02 


23 


37 


| 2.79 


! 1.31 




15 


06 | 


14 


4 6 


94 


| 3.04 


I 1.97 


1 8 


33 


32 | 


03 


12 


100 


| 3.20 


I 2.08 




15 


56 | 


15 


36 


134 


| 3.33 


1 2. 16 


1 9 


34 


19 | 


03 


59 


133 


I 3.46 


1 2.25 




16 


42 | 


16 


22 


113 


I 3.59 


! 2.28 


. ... _ 












j 


j- 


20 


35 


05 | 


04 


45 


113 


| 3. 50 


| 2.23 




17 


26 | 


17 


06 


139 


| 3.48 


1 2. 26 



and 



2 1 
22 
23 



35 47 
13 03 

36 28 

18 49 
3 7 09 

19 33 



at 

•- — 4 



the shelfbreak Hiah 



05 27 
17 43 

06 03 
13 29 
06 49 
19 13 



I 



135 

132 

97 

93 

33 

76 




maximum density on the eastern half of the ground coverage, 
as shown by a plot of the spatial-temporal distribution 
(Fig. 4.7) with respect to the distance from the shelfbrsak 
(y-coordinai e defined in 1-C and Fig. 1.1). 

1 * Analys is of the Mean 



A smoothed transect was established by averaging 



along the x-axis with a " y" step increment of 
4.8) . The isotherms generally sloped down from 
region to the open sea (on the average 10M over 
marked mixing was apparent in the region of the 
(-10<Y<10) with shoaling of the isothermal layer 



2 0KM (Fig. 

the shelf 
100KM). A 
s hel f br ea k 
s above the 



thermocline (14°C) and sinking below it. Dynamic heights 
relative to the smoothed transact yield a mean surface gecs- 
trophic flow relative to 150i1 of the order of . 2M/S and 
about 20KM wide, towards the northwest at the shelfbrsak. 
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TABLE 71 

SURFACE TIDAL CURRENTS AT H7°24'N 006005'B 
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The geometric approximation of the shelf by a 
straight lire in the averaging process, and probable varia- 
tions in time and space of this mixing region, reduce the 
more dramatic features that can be seen in oarticular non- 
averaged transects. The mean vertical -emperature profile 
is characterized by a mixed layer at 130c limited by a sharp 
ther mocline (on the average 0.2°C/M| . On the shelf, a bottom 
mixed layer of 12°C was present. 

Taking an arbitrary "level of no motion" at 150M, 
dynamic heights were computed for each X3T soundings and a 
horizontal smoothed contouring was made (Fig. 4.14). The 
surface gecstrophio speeds corresponding to this analysis 
are cyclonic around a stretched ellipse with a major axis 
coinciding with the continental slope. Northwest currents 
or. rhs uroar oart cf tha slooa ara of . 5 H / S at fhsir 
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maximum, vs. a returning southeast flow of . 23/S over the 
rise . 

A vertical modal analysis (Cf. 71-3) , yields baro- 
clinic internal Rossby radii of deformation of 4.5KM on the 
shelf and 25K.1 in the deep ocean. ( 3 f VI-B) 

2. Analy st s of the Variability 

The spatial variability of the thermocline depth was 
computed using a subset of the total X3T data, which was 
selected to avoid unrepresentative (sparce and dispersed) 
values. The variance of the depth of a given isotherm was 
calculated in each rectangular bin of IQFM width. Standard 
deviations of the 1 3°C isotherm depth from the above means 
were computed (Fig. 4.9) . The strong increase in the values 
at the position of the idealized shelf break could either 
mean a significant increase in the time variability 
(internal waves), cr express the fact that the mixing is 
taking place over a very short distance compared with fh a 
width ( 1 0 KK) of the bins used for averaging. 

3 . S el ec t °d Rep res s r.t a tl ve Transects 

Seme specific transects were then considered to give 
a quasi- instantaneous spatial image of this mixing process. 
One interesting subset of 3T's was acquired by a transiting 
ship sailing towards Ouessart (Osh ant Island) and whose 
track crossed the s’nelfbreak at 47O30N, 95°W on August 23. A 
sampling rate of one per hour (i.e. 30K.d spacing) enables 

drawing a quasi-synoptic transect parallel to the SAR swath 
but 60KM west of if (Fig. 4.10). Other transects were made 
by selecting 3T's close in time and conveniently spaced in 
the vicinity of a perpendicular to the isobaths. A se~ of 
three transects is depicted on a geographic grid (Fig. 4.11) 
with other non- synoptic transects reproduced from various 
sources (Fia. 4.12) . The scaling adopted is identical for 
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all and bas=d on the "alor. g-^rack" distance from the shelf- 
break. Common characteristics appear to be a narrow pelagic 
front anchored on the sea-floor at the upper part of the 
shelfbreak (ca . 180M depth} , with one or two upwelling zones 

adjacent to this narrow intermediate mixed layer. The ther- 
mocline is stronger .on the shelf than in deep waters (Fig. 
4.13). In chapter VI I , an idealized 3-layered ocean was 
extrapolated from these selected transects and the corre- 
sponding geostrophic currents were compared to the observed 
circulation . 

E. TEMPORAL THERMAL ANAL? S IS 

Daily averages of the isothermal depths (Fig. 4.15) show 
a significant heating in the upper layers on 20 Auaust and a 
general broadening of the t her modi ne between 17 and 24 
August. Assuming a propagation offshore from the shelf of 
tidally generated waves, diagrams have been established 

displaying for a given isotherm, variations of its depth as 
a function of time and space (Fig. 4. 16) . These variations 
in depth were obtained by subtracting the mean along the 
x-directicn to eliminate any permanent trend and enable a 
contouring of relative values. Tne spatial axis is th® 
along-track (y-coordinate) distance from the shelfbreak; 
time axis in hours has osen "folded ov®r" two semidiurnal 
(12. 4H) periods. Hence all time coordinates are *aken 

within 24 hours by subtraction (or addition) to the actual 
time of observation of an adequate multiple of the semidi- 
urnal period. 

Mo contouring algorithm could be applied to such sparse 
and greatly contaminated data (by various phenomena; vari- 
able wind stress, radiation heat exchange, even high 
frequency nonlinear internal waves). A hand-contouring was 
attempted under the following or i n ci pies : one contour line 
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was drawn fcr a zero relative depth, another line followed 
the maxima and the last one followed the minima. This was 
done for three temperatures (12 , 1 4 , 15 °C) . Each of ■ u hcs- 
diagrams exhibits a constant pattern obliquely oriented on 
the axis, and thus corresponding to a propagation of the 
variations. The slope of those patterns (disrates of propa- 
gation during a given interval of time} was found to be 
approximately 1.1M/S (equivalent to D.9M/S perpendicular to 
the shelfbrsak). A modal analysis described in VI-6 showed 
that high frequency internal waves of period 90iMIN would 
have this same wave speed. Furthermore the phase of -his 
propagation implied a maximum deepening event an ,, Y = 0'’ at 
about 1900Z, 20 Aua. Conversely, an the time of nhe SAR 

image the crests of these propagating oscillations were 
found to be at Y=10 and Y=60. Maximum values of relative 
depth excursions show general increasing values with depth: 
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At deeper depths, the MEADS mooring temperature data 
allows estimation of the vertical isotherm movement by 
relating semidiurnal temperature fluctuations at a given 
depth to the cl imato logical temperature profiles. Data from 
the current meter pressure gauge showed that vertical oscil- 
lations of the array itself were negligible. 

At 1500M, temperature variations of .3°C in a mean gradient 
of 9 10* 3 °C/M yield vertical amplitudes of about 2 0 ‘1 . 

At 3000M, temperature variations of . 15°C in a mean gradient 
of 2 10'* °C/M yield vertical amplitudes of more than 40M. 

At 3500M, temperature variations are at the limit of 
discrimination of the thermal gauge but reveal vertical 
movements of at least 20 M. 
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F. SUHHARY 



SEASAT flew over th® shelfbreak of the Say of Biscay at 
2142Z (after sunset) cn 20 August 1973. A high pressure cell 
had dominated the area since early in the morning, brinoina 
a southward wind of 651/S. The cloud cover gradually 
decreased, and the air temperature decreased to 17°c. 

A well established swell was entering th® Bay from the 
northwest, nearly tangential *o the isobaths in the vicinity 
of the shelfbreak where a significant increase in wave 
heights was observed from the SEAS AT altimeter. 

The tide was a spring tide at high level plus 5 hours at 
the shelfbrealc. The ebDinc surface current on the shelf 
side was at its maximum toward 240°r at 0.8551/S. 

The ocean was strongly stratified with sea surface temp- 
eratures of 18°C, a thermocline gradient of 9.3°C/51 at an 
average dep^h of 30M, and a bottom iidally mixed layer on 
the shelf of 12°C. The first mode internal Rossbv radii of 
deformation were of 4.50 or. the saelf and of 25 0 in deep 
water . 

A spatial thermal analysis showed a mariced front at th® 
shelf break with intense mixing and upwelling. There was a 
westward surface geostrophic get at the snelfbreak. 

A temporaL thermal analysis gave some evidence of an 
internal semi-diurnal "disturbance" propagating offshore 
from the break with a speed of D. 9 51/S and with the last 
crest in the vicinity of "Y=10" an the time of the 3AR 
image. This propagation sp-ed being sigr.if ica r t ly smaller 
than the computed phase speed of the internal first modes 
tides (3.8 and 1.751/S), either some higher modes are presen' 1 -, 
or high frequency waves are generated by the tide but even- 
tually decoupled from it and trapped in the vicinity of the 
thermocline (where the largest Bru.nt-Vaisala frequencies can 
be f o un d ) . 
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Time. (Aug-2.o) 1*4 



7. SAR IMAGE ANALISIS 



A. ORBIT PARAMETERS AND GEOGRAPHIC REFERENCES 

Tha auxiliary grid introduced previously (1,4) will be 
usad to cor. van lent ly locate features with respect, to tha SAR 
image and to the average shelfbreak contour. Tha geographic 
fit of the SAR image is estimated to be accurate to a -1- least 
+/- 5KH in tha vicinity of the shelf-break. 

At 2142Z Aug 20, SEASAT was on a descending rode over 
France with the SAR antenna oriantad to tha northwest. 
Consequently, a positive (northward) shear of eastward 
current would appaar on the image as artificially compressed 
with enhanced Bragg-scatter (Cf Fig. 2.3) In oth<=r words, a 
current directed toward the northwest would appear on th=> 
SAR image with a bright northeastern boundary and a dark 
southwestern one. 



TABLE 711 

SEASAT ORBIT 785 PARAMETERS 



Orbit 7S5 on 


August 20 1973 


(Julian Da 


y 232) 


T 

1 

— 1 


Ascer.dina node at 


[ 2103. 55Z 






East Longitude 


| 168.43 02 


( deci mai 


d = g r = e s) 




Descending node at 


I 2154. 14Z 




East Longitude 


I 335.39 58 


(decimal 


d ears e s) 




Argument of Perigee 


| w=140 




Recent rici tv 


1 e=0 .0016 
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HEURISTIC A HALTSIS OF THE SAR 1.1 A GE 



The patterns visible in the optically processed SAR 
image may be grouped in four distinct categories according 
to their geographic location and to their geometric charac- 
teristics. They are discussed separately under the titles: 
"inner shelf", "outer shelf", "generation area", and "propa- 
gation area". The image discussed being a positive print, 
the lowest backscatter radar response will be referred to as 
dark areas and, conversely, the enhanced signal returns will 
be referred to as bright or light areas. 

A correlation between image char act eristics and bathyme- 
tric features is made and must take into account two inaccu- 
racies. First the error in mapping the SAR image is 
estimated as 5KM at most. SecondLv the bathymetric chart 
used here (Fig. 5.5) was constructed by interpolating sound- 
ings of different sources, with gaps of ca . 10KM in seme 
places of the continental slope. Hence, the relative posi- 
tions and shapes of the canyons drawn are subject to 



caut ion 


within the 


same order of 


magnitude 


as the SAR 


feature 


locati on s. 


Nevertheless, 


in the case 


of several 


assumed 


cor rela tion s 


between well-def 


in ed imaae 


and bathyme- 


trie f 


eatures, the 


confidence may 


be increa 


s sd by the 



perfect relative fit of each of the points (Fig. 5.4). 

1 . Inn er Shelf 

(Y>40) - In this area, manifestations of internal 
activity are of low contrast and lack any obvious organisa- 
tion. Separated curved lines can nevertheless be identified 
with an average radius of curvature of 12 KN. Patterns of 
similar shapes can be associated at spacings of 28 0 (Fig. 
5.2) , but their relation to any kind of periodic phenomena 
is net proven. They could be due to mesoscale atmospheric 
variability, or linear slicks associated with a localized 
convergence (If para n-c-2) . 
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The lowest limits of wind speeds sufficier. 4 - for a 
radar return are attained in the northern par- of this area, 

l 

particularly south of the "Chausses of Sain" and leeward of 
the Britanny Toast where mesoscale land-sea breeze may well 
have been predominant. (A ship sailing 100°T can clearly be 
seen with its wake . virtually displaced according to the 
Doppler effect described in II-B-li. 

2 . Outer Shelf 

(10<Y<40)- Internal wave group patterns are apparen 1 ’ 
in this region. They appear as arcs of circles, cor.cav= to 
the west with a variable radius of curvature of minimum 
value 15KM. They are limited to the north by a caustic in 
which they merge; this line is brighter, i.e., of greater 
backscatter. They are broken tc the south by another line 
of brighter characteristics still. These patterns seem to 
exist over the entire width of me 5AR image, but they 
appear with more contrast on the western par- of if. Som = 

periodicity, of about 29KM lengtn on the x-axis, can be 
detected . 

The wavelength of these patterns is variable in two 
directions. (1) toward the so-called caustic to the north, 
the wavelengths decrease to zero, starting at a maximum 
value on the southern boundary. (2) toward the southeast 
(increasing x-directior.) , which is their apparent direction 
of propagation, the average wavelengths, measured at mid- 
distance between the two boundaries, is of the order of 1KM 
on the western side, increasing to 2K.1 in th a middle of the 
swat h . 

The most conspicuous wavecrests are located at 
x=-25, where a high level contrast particularly enhances 
those curved lines. Bach of these arcs appears as two 

narrow contiguous stripes: a dark (low backscatter) one on 

•he concavity and a bright (high backscatter) one or. th = 
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convexity. Hence, over one wavelength on a dir-c~ior. 
perpendicular to the crests and following the assumed direc- 
tion of propagation (from the inferior to the exterior of 
the arcs), the backscatter varies approximately in the 
following way: 70 % of the length at a uniform level, 20% at 

low backscatter level and 1 0 % at a nigh level on the leading 
edge. According to the hydr ©dynamical theory of SAR imaging, 
this result means that the internal waves responsible for 
these patterns are long-crested with deep narrow troughs 
(Fig . 2.4). 

The southern boundary is a general line brighter 
than the ncrtnern caustic particularly on the western part 
of the image, where it is conspicuous both by its brighter 
characteristics and by the obvious discontinuity it creates 
in the general patterns. This boundary is located approxi- 
mately (+/- 5KM) over the 180M isobath, thus corresponding 
to the upper break in the continental slope. 

3. Continental Slo pe o~ Ge ne r at ion Area 

{ 0< Y<1 0 ) — This area covers the upper part of the 
slope between the 150 and 1000M isobaths. Once again the 
mesoscale features are concentrated o.n the western side of 
the image. The average backscatter is weaker (darker 
shading) than on both adjacent areas (shelf and deep ocean). 
The northern limit is the band of lighter characteristics 
described above. 

The limit to the south is (conversely) a discontin- 
uous band of very low backscatter about 15001 wide, which is 
clearly visible from x= -6 0 to x=-20 and from x = -5 to x= 1 5 . 
The straightness of this line (vs. the other patterns) and 
the fact that it is not bordered by anv alternating clearer 
stripes, seem to indicate that it is due, not to a propa- 
gating phenomenon, but to a stationary (if not steady) 
internal feature. These could be symptoms of a surface 
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divergence according to the kins metical theory of SAP 
imaging. 

In this area, well-orga.ii ssd coherent circular 
patterns suggest cylindrical spreading of internal waves 
from several distinct point sources. Assuming that the 
direction of propagation of a given wavepacket can be esri- 



mated both by 


the curvature of 


the crests 


and by th- 


tive vaveler.gt 


h dispersion sense 


, the fo 


llowing mai 


souces can be 
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for further 
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depth . 








Point 


P2 is at the head 


of a narr 


ow and obli- 


deep, unnamed 


ca nyon . 






Point 


?3 is at the head 


of Ouilvi 


nec Canyon : 



dept h . 

Point P4 is at the head of 31avst Canyon at 300M 

depth . 

It is observed that, at the Oder canyon there is no 
evident corresponding internal wave pattern; a plausible 
reason for this could be that this canyon thounh steep, does 
net go up to the shelfbreak i*self (Fig. 5.5) . 
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At PI the leading wave cf the nearest group is 15KM 
from the center. Wavelengths from the leading *rough 

decrease from 2 K M tc 0.7KM with a range of +/- 50*. Som = 
interference can be seen with a second circular pattern of 
poorer definition centered around ?’1 ( X = - 4 4 , Y = - 1 0) . 

At P2, two wave groups propagate in two distinct 
directions. The first group propagates toward 130°T along 
the X-axis with its leading wave at 10KM iron the center. 
The second group propagates with a 33° inclination offshore 
(160°T), its leading wave being at 18KM from the center and 
having the greater wavelength in the group as above. 

At P3, no obvious group can be identified ba 4 the 
pattern appears as a superposition of faint ellipses with 
their major axis oriented 130°I along the X-axis. 

In the vicinity of P4, patterns are still more 
confused, but a guasi-circular group extends around it at a 
distance of 17 KM. 

4. Continental R ise or Prop aaa 1 1 on Area 

(-90<T<0) The surface patterns hers are mors blurred 
than in the three areas described above. Features are on a 
larger scale and do not present the same east-west vari- 
ability. The brightness is uniform over the whole area. The 
general geometric organization is that of arcs of circles 
roughly centered on the point sources defined above, or more 
generally, on the 1Q00M isobath. Some alternating black and 
white bands may be identified, but are larger and do not 
contrast sharply on a grey backdrop (in contrast to the 
so-called generation area). The distribution of those 
circular bands is not always clear, but some particular 
patterns with similar shapes can be found at a spacing of 
ca. 23 KM. 
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ANALYSIS DP THB SWELL 



The remotely generated swell is long and well organized, 
thus giving generally well identified patterns on the SAR 
image. Some apparent geographic variability in direction and 
crest length can be observed and could be indicative either 
of variations in the nonlinear imaging mechanism, or of real 
deformations in the swell field by interaction with surface 
currents or mesoscale winds. The wavelength of the swell is 
237M (an 47°N, 00 4 ° . 3 0 1 5J) +/- 30 M. This corresponds to a 

period of 12.3 seconds, a phase speed of 19.2 M/S, and a 
group speed of 9.6 M/S. The direction of propagation is 



well defined 


and 


easily measurable 


d n 


the shelf and 


offshore 


Spain on the 


s 0 u t 


hern part of the 


UJ 

UJ 


image. In the 


vicinity 


of the Armor! 


.can 


shelfbreak, the 


pit 


terns are nore 


con fu sed 



and mix with the infernal wave field. On the shelf, the 
general orientation is 115°T (at 47<JN, 004OW) . In the open 
ocean (46°N , 006° W) , it is 135°T. South of the shelfbrsak 
if seems to drastically bend towards the south as if 
refracted by an opposing current. This particular feature 
will be discussed further in sectioi VT-E. 

D. SOHMARY 

Or. the shelf, and, particularly in the first 30 KM 
inshore of -‘•he she 1 f bre ah, the general orientation of the 
surface patterns suggest an internal motion propagating 
along the isobaths and toward the east. A strong and local- 
ized divergence follows '-he shelfbreak at 180M of depth, 
clearly limiting to the north a generation area on the slope 
where several point- sour css can be identified at -he heads 
of sharp canyons. Internal waves are generated in packets 
which are clearly visible at short distances from their 
generation points. Wavelengths of the internal waves, as 
well as the spacing between packets, are highlv variable, 
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VI. CORRELATIVE INTERPRETATION OF THE SAR IMAGE 



The entire area cf the continental slope imaaed by SAR 
is dominated by surface manifestations of shore infernal 
waves. Thsss internal waves appear in. groups which are prob- 
ably generated at the shelfbreak in association with a 
semidiurnal internal tide, but eventually decouple from it 
as both perturbations propagate offshore. Two different 
fields of internal waves have to be considered; (1)the 
internal tidal wave with a semidiurnal period of 12. h hours, 
and (2) the higher frequency waves generated in packers as 
forced by the internal tide interacting with shelf bathy- 
metry, or by dissipating energy at higher wavenumber by 
breaking processes. The c haracteri s tics of these internal 
waves are derived from their surface manifestations (Cf 
chapter V). Comparing the results with a synoptic conven- 
tional analysis (Cf chapter IV) of thermal structure and 
Eulerian current, measurements, gives a notion of the 
validity cf an operational approach where a SAP. image could 
be the basis for hindcasting the actual internal motions. 



A. TWO LAYER MODEL ON THE SHELF 

The sharpness of ‘■he thermocline on the continental 
shelf allows the assumption, in a first step, of a two-layer 
ocean and use of the corresponding dispersion relation to 
compute phase speeds for given observed wavelengths. The 
dispersion equation is: 
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Or. the outer shelf, for a giver, internal wave crest, the 
observed increase of wavelength from the "caustics" {Cf 
V-B-2) to the leading middle point averages 25ft. Working by 
logarithmic differentiation on the simplified form of the 
dispersion equation 

c i = _ap_ Wh-M 
j e h / 

gives an approximate relation between this wavelength rela- 
tive variations and the increase in mixed layer depth: 

l - 1 ^ ^ ^ ^ 

C A \\ V\ — V\ 

and 

Ui , i iL ^ (h-M . 

A U - ZV\ 

According to the observed variations and to the mean 
thermocline depth derived from X3T data, -he thermccline 
depth is estimated to vary on the outer shelf from 25‘1 on 
the boundaries to 55M at mid-distance from each of *h c o. 

On the continental slope, a linear interfacial wave of 
wavelength 1KM yields a phase speed of 0.48'1/S whereas for 
700M the phase speed decreases to 3. 47. M/S. Applyir.a these 

results to the distance separating the wave packers gener- 
ated on the shelf , and assuming one generation points are 
those bathymetric features described above (Cf. III-B) , the 
travel time would be 9 hours for the leading waves, thus 
implying a generation at low tide with maximum ebbing 
surface currents (Cf Table 4.6). 

Conversely, and if the generation time was known with 
better precision relative to the tide phase, observation of 
the position of similar groups together with density 
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gradieno data could give ar. estinace of oh? *■ hermo cl in ? 
depth . 

This two-layer approximation would hold if the chermc- 
cline was sharp enough or equivalently if, ao its level, the 
Brunt-Vaisala frequency profile exhibited a strong peak, 
thus trapping the energy of high frequency internal waves in 
an horizontal channel. But indeed, the actual ocean has not 
such a simplified structure , particularly in the "genera- 
tion area", where dramatic mixing processes are talcing 
place. Furthermore, at low frequencies such as semidiurnal, 
the vertical movements extend to tie whole water column and 
cannot be described without a vertical modal analysis. 

B. NORMAL MOOES OF THE BAROCLINIC INTERNAL TIDE 

Internal baroclinic modes of the ocean were computed 
from the density profiles extracted from observed temoera- 
ture soundings and climatological salinity data (see Table 
VIII) . The Brunt-Vaisala frequency is the free oscillation 
frequency for a given density gradient : 

n 2 = - J_( . 

e ' 42 c 1 > 

Density gradients were calculated by o he equation of soate 
from the observed temperature profiles sampled by X3T data, 
and assuming a constant vertical salinity profile of 35. 6^0 
in the deeD water and of 35.u%c on ohe continental slope 
from climatological data. Bruno- v a isala frequenci-s wer® 
then computed in the first 300M and eventually blended 00 a 
deep mean climatological profile found in Gould * 1973]. 

The eigenfunctions solving the Sturm Liouville problem 
of veroical modes give the vertical dependence of the 
velocity fizLslds Ho zr ^.ny "zr^crus'i^y oz. z.^.2. 
oscillation. An infinite number of soluoions exiso yielding 
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vertical distributions characterized by one or more phase 
reversals in horizontal speeds witn depths (Fig 6.1 and 
6.2) . For eacn of these modes internal wavelengths and phase 
speeds can be found as a function of the fencing frequency. 

Computations were accomplished using a numerical solu- 
tion [Mooers 1971] to the eigenvalue problem. Starting with 
the boundary condition a 4 - the surface, successive correc- 
tions cf the eigenvalue are made at each iteration to 
satisfy the bottom boundary condition. 

The solutions are defined «/ithir. a multiplicative 
factor. An estimate of the vertical distribution of ampli- 
tudes values is then necessary to resolve this incertitude. 
A simple best fit of the amplitude values from available 



data was made to the first baroclinic mode. Vertical ampli- 
tudes at the semidiurnal frequency were extracted from the 
current meter data (Cf I V - C ) : temoerature fluctuations as 

related to the mean vertical thermal gradient at the corre- 
sponding deDth gave an estimate of the vertical excursions 
of the isotherms (as shown in IV-Ei . (The variations of the 
current meter pressure gauge showed that vertical oscilla- 
tions of th e array itself were negligible) 

The second mode vertical profile is coherer.* with those 
observations and is consequently assumed to be the mean 
component of internal energy at the semi-diurnal frequency 
(Fig. 6.3). Maximum vertical amplitudes are inferred to be 
a~ 2000M and of the order of 50tM, a level of no vertical 
movement (medal node) existina at 1333M. 

The analysis of the 31 soundings exhibits an offshore 
spacing of perturbations of 45KM corresponding to phase 
speeds of 1 M/3 for a semidiurnal generation and wavelengths 
of 7.2KM in the open sea (90MIN period! from the first mode 
dispersion analysis. 
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I?, conclusion 



9 



~ .C « i 

it j o 



, the inte rnal tils, at least in. 
mode, does not appear an the S A3 image, neither in the 
spatial analysis of the the rmostruct urs . Nevertheless, peri- 
odic disturbances could be generated at a semi-diurnal rate, 

consisting of high frequency waves which are discussed now. 



TABLE VIII 

INTERNAL WAVELENGTHS 



WAVELENGTH (KM) 





PERIODS | MO DEO 
(Hours) | tropic 


M3 DEI | 
vs. 


| MODE 2 

3a roc J 
! 


MDDE3 iMODE4 
.inic modes 
1 


SHALLOW WATER 
1 SOM 


12.4 j 2 5 70 . I 
6.0 | 8 88 . 

1.2 1 166. | 

0.4 j 55. 


27.3 0 
9.30 
1.73 
0.53 

4 3 To ~ 
14.9 
2. 7 
0. 8 


10.00 | 5.60j 3.50 

3.50 | 1.90| 1.20 

0.63 1 0.351 0.22 

0.19 J 0.121 0.08 


MIXING AREA 
2 20M 


12.4 13110. 

6.0 |1076. 

1.2 1201. 
0.4 | 67. 


h 18 . 00 |~ 12.07) 3.98 
6.22 | 4.17| 3 .1 0 

1.10 1 0.76 1 0.57 

0.30 1 0.21 j 0.16 


DEEP WATER 

3500M 

J 


12.4 114080. 

6.0 |4870. 

1.2 1 382. 

1 1 I 


"T7o7o~ 

53. 4 
2. 3 

1 A 


^76 . 0 144.00 {29.00 

25.8 |15.30 | 9.70 
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HIGH F3EQ0ENCY INTERNAL WAVES 



The wavelengths in ths internal wave packets are of the 
order of 3KM, yielding frequencies of the order of 70 MIN 
which are consistent with the internal wave climatology (Cf 
III-F). The spacing of the high frequency wav* packets as 
seen in the propagation area of the SAR image is approxi- 



mately 23KM. 

On the continental shelf, the periodic spacing in ths 
patterns is about 28KK, which is olose to the semidiurnal 
internal wavelength computed for shallow water. 

1 • E stima tion of the Nixed Laver Depth 

It is assumed that the climatological density struc- 
ture of the Bay of Biscay is known, and that frequent XBT 
surveys as well as remote sensing of the sea surface temper- 
ature, allow estimation at any time of the mean density -jump 
a- the thermocline. Starting wifi only this hypothetical 
element and making abstraction of the actual in-situ data 
available, it is attempted to compute the thermocline depth 
as a function of the geometric characteristics of internal 
waves as derived from the SAR image. 

Assuming a climatological density jump of 1.5x10 3 
at the thermocline, and observing separation between wave 



packets of 23KM in the op=n sea area of the SAR image yields 
an estimate of the mixed layer depth of 20 M. A variation of 
10 % in the assumed density jump would result in 10 v - differ- 
ences in this estimate, whereas changes of 10 % in the inter- 
packet distances would result in 20^ changes in this derived 
mixed layer depth. Consequently, if the absolute estimate 
of the thermocline depth is uncertain, comparison between 
interpacket distances would be a very sensitive way of esti- 
mating time variations. 
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• Est i m at ion of th r I r.t e r nal Wave s Am pi it ud a 

Amplitudes of the nonlinear short internal waves 
were computed by a method described and applied cn a similar 
SAR image by Apel and Gonzalez [1993]. The firs- assumption 
made was that in the "generation area" (continental slops), 
the internal waves visible on the SAR image could reasonably 
be considered as solitons. By definition, a soliton is a 
perturbation conserving its shape and general characteris- 
tics (except phase) , even after collision wi'-h another one. 
This characteristic is apparent in the vicinity of point Pi 
for example (Fig. 5.3), where several wave packets are 
superimposed without any notable mutual interference. 

This conservative nature of solitons is linked to 
the fact that the linear phase speed dispersion is counter- 
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wa v e 1 an g 
ir.t ernal 
tion to 

was defi 
linear t 
1933], 
elliptic 
argument 
rude A: 



th, and a phase speed limit "C " for infinitely long 
waves in a non rotating ocean (steady-state sol ti- 
the Star m-L i ouvi lie problem) were computed. 



In one 


given 


packet 
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the 


wavelengf 


h st re 


fch, (»f'') , 


ted as 
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ratio o 


-C 


the 


observed 


wa vel 


ength to the 


heoret i 


cal o 


ne . F r 


om 


the 


theory 


C Apel 


and Gonzales 



this ratio is equal tc "IMm 2 )", twice the complete 
integral of the first kind with argument m 2 ; the 
m 2 beina itself a function of the sciitcn's ampli- 



c °^- -A- 
“ C ? ^ 

For any wavelength if is then possible to extract 
the nonlinearity due to a given amplitude. dr, conversely, 
observing a nonlinear dispersion in the patterns of the 
internal waves (by measuring the variation of wavelength in 
a given packet) , if is possible to derive an estimated 
amplitude of the waves. 

was taker, to be of the order of 5 :c 1 0 ~ 2 m l , ("oC’* ” = 20h) . 
”C 0 " was found to be 0.55.d/S. 

For each observed wavelength in a given packet, the corre- 
sponding phase speed can be computed, and the wavelength 
stretch can be measured. A value for "m 2 " is derived from 
which wave amplitude is calculated as 



(\ - c ° ~ c P 3 m“ 
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ng waves were found, 
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D. INTERNAL SATE INTERACTION WITH 5 H2LFB2SAK STUCTOSE 



Two different kinds of geographic changes in the 
internal wavs group characteristics are apparent (Cf 7-5-3). 
The general shortening in wavelength will be assumed here to 
be either a Doppler shift due to a current shear or a phase 
speed compression due to a horizontal density gradient (e.g. 
as in a front or a mixing area) . Refraction of the crests 
resulting in a deformation from their original cylindrical 
spreading as in groups northeast of ?1 and ?2 , can, on *.h« 
ether hand, be more surely associated with a strong near 
surface current. Estimates of the direction of a shear 
surface currant as given by rh aracteristic refraction 
patterns are indicated on figure 5.4 They give some indica- 
tion of a general southwest ward alongshore flow no the south 
of the generation area (on the continental slooa) , opposite 
to the more localized currents at the shelf break, and to 
those south of the slope (Fig. 5.4). 



E. SWELL REFRACTION AND S3 ELF3REAK CURRENT 



A subject 
attempted. No 
but in most pa 
such that visa 
precision) cf 
“sting to try 
interaction 
refraction is 
probable. T 
[19 8 1], has b 
by Neado ws [ 19 
caticns to thi 



ive analysis of tns swell refraction was 
Fourier Transforms were mads of the SAR image, 
rts of it, the clarity of the swell field is 
a 1 estimation of its direction and (with less 
its length can be made. It was thought inter- 
ne apply recent theories on wave/current 



o this particular case, where some severe 
obvious and the oresenre of an interacting jet 
he theory initially developed by Phillips 
sen applied to a SAR image of the Gulf Stream 
331; it was transposed with drastic simplifi- 
s ca s e . 
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In the same way that the Snell's law (loi de Descar~.es), 
describes the refraction of a rav encountering a era die r.t of 
propagation in a medium, the following equation describes 
the refraction of swell of a given period, propagating with 
a phase velocity "c" through an interface between owe sea 



surface areas 


with surface 


currents of " u 0 " 


and " U. , " ; 




an d " o( , " ar e 


the angle of 


incidence and of 


refract icn 


o f 



swell "ray" (Fig. 5.5). 



Cos (X, = 



CoS X. 






-\z 



Assuming that the cause cf the observed refraction is a 
shear current directed along the shelf break (cr at least 
parallel to the isobaths), the angle of incidence and the 
angle of refraction (just after the first visible break in 
the wave crests) were measured. Tie current speed (equiva- 
lently the get velocity) was estimated to be 0.5J1/S. The 
accuracy obtained cannot be more than 50^, but the basic 
observation seems reliable enough , at least, to suggest the 



oresence of a 



jet and to localize 



he 



oosition or 



the 



maximum soeed at the continental break. 
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71 I. TENTATIVE EXTRAPOLATION £0 GENERAL OCEANIC FEATURES 



A. I STERNAL WA7E GENERATION 



"In the ooasta 1 reg ion internal 
continental slope are standing waves 
modal analysis relative to the iensit 
the tine of the SAR image shows wave 
27 KM for the semidiurnal first mods 
shelf. This wavelength is comparab 
patterns visible on the SAR image no 
If can therefore be assumed that t.os 
to the interaction of another interna 
cally modulated in phase with ridg 
internal fids. Obviously, the as 
standing internal fide would only be 
of the shelf sea acted as oerfect ref 
Wavelengths of the first intern 
deep water are 17 0 KM for mode one, 
44KM for mode three. No direct 
length scales is visible on the SAR i 
high frequency internal wave packets 
st rat ion cf the high energy loss enco 
fide at the irregularities in the she 
heads of the narrow canyons, osri 
column are estimated (for this speci 
tide) to be of the order of 15M wnir 
waves forced by the ebbing tidal curt 
packets generated by this mechanism 
canyons extending toward the top of 
1000 a isobath and occur in bursts, 
an isctro cicallv due to the oresence o 



tides generated on a 
'• [Rattray 1959 ]. The 
y structure observed at 
lengths cf the crier of 
in* ‘Sin si or. r h e 

la with th? soacira cf 
rth of -he shelfbreak. 
se patterns correspond 
1 oerfurbafion oeriodi- 
=5 and troughs in an 
sumption of a perfect 
true if fc he boundaries 
lectors. 

al barer liric modes in 
75KM for mods two and 
manifestation at these 
mage. Nevertheless, the 
give an indirect demon- 
untsred bv the internal 
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are localized at the 
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They eventually develop 
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the north. The general direction of propagation is toward 
the south-west perpendicular to the continental slope. A. 
thorough discussion of the possible generation mechanisms is 
given by Pingree [1981]. The specific influence of narrow 
canyons has been recently enlightened by a laboratory exoer- 
iment [Baines 1 983 ].. 

B. TIDAL MIXING ON THE CONTINENTAL SLOPE 

Tidal nixing on the continental slope is demonstrated by 
the divergence in isotherms in the quasi-synoptic X3T anal- 
ysis, (Fig. d.10 and u.13) . and is in good agreement with 
other measurements made in summer and fall at the same loca- 
tion. In Pingree (1981), the various hypothesis related to 
this mixing are developed. 

C. SHELFBREAK FRONTS AND JETS 

A new insight on the sheifbreak circulation is given by 
the SAR image. Parallel bands of divergence-convergence 
which compress internal waves, can be associated with the 
thermal structure inferred from BT data and illustrated by 
the following simplified scheme of circulation (Fig. 7.1). 
Over the shelf break, the dynamic heights tend to be lower 
than on either side as a consequence of the tidal mixing. 
Therefore, a strong geostrophic flow extends to the therno- 
cli.te, follows the isobaths, and is confined to the shallow 
part of the continental slope. Tie flow is consequently 
toward the northwest on the upper mart of the slope and 
toward the southeast in the off snore half. The idealised 
transect (Fig 7.1) , accounts for the localised and strong 
northwest flow and explains, by stretching of the thermo- 
clin e o 5s shears, ^hs icwsr vaiuss for "h.9 co u nt * 

cu rr<9nt. I* is nevemhsless o b 7 i 2 as “ha* 3. n v ’<ind of 
analytical model oased on geostrop h y (steady-state) and 
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w he r. 



alongshore homogeneity cannot be fully appropriat e 
applied to the Celtic Sea s he 1 threat . In this area currents 
are dominant and limited geographically by the coast of 
Britanr.y, bathymetric irregularities are essential in the 
internal generation processes, and nonlinear effects between 
the steady-state circulation and the tidal movements may be 
expected. 



The most informative indication 
strong alongshore let are the swe 
sudden increase in wave height (IV- 3, 
cidence of these phenomena give reli 
which may be related to the pole 
current discussed by Swallow [1977], 



of fh9 Dr9S9n.C9 cf 3 
1L z 9 f z 9 cf i z n f. oh 9 



7-C, 71-3) . 

3 b i 9 9V3333C9 

:a ^ nr r 



7 h9 coin- 
of a jet 
boundary 



106 



r 



-=> VJoQxH 




Aays3Au 



The barotrooic semidiurn al tide Drooaoatss onshore, 

generates a barcclinic internal tide,' 

eventually reflected by the steeo continental sic oe . 

Another barociinic tide is generated on the shelf, 
and reflected both at the coast and the shslfbrsax. 

At the canvons, traooinc of internal energy induces 
larce aaolitudss vertical aoveaents, ceneratinc 
high frequency waves, propagating on the t her accline 





Surface a vn a. ale heiohte and t^ostroohic flow. 



•figure 7.1 3C32KA2IC IY?0 THESIS OF SHILI3S2AK CIRCULATION. 
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VIII. APPLICATION TO ONDERffATER ACOUSTICS 

Of the host of variables affecting sound pcopagat ion in 
the ocean, the temperature profile has the strongest influ- 
ence. Not only the vertical gradient profiles, out also the 
horizontal gradient fields, determine the widely varyir.u 
shapes of the insoni f leaf ion volumes. 

Hence temporal and spatial temperature variability 

should be taken into consideration when forecasting any 
acoustic quantity, from one single climatic or observed 

temperature profile. Two approaches are possible in ■ his 
matter: (l)a deterministic approach would be possible only 

in the cases where the water-mass's vertical movements are 
known as functions of time and geographical locations, well 
enough to allow the forecast of effective variations of the 
acoustic transmission loss with time, location and orienta- 
tion of th* sound • beam, and (2) a pro ba bid is t ic approach 

would attempt to model “he oscillations without any phase- 
lock and would only yield estimates of the acoustic propaga- 
tion variability. 

Obviously, in the common case where the geophysical 
variables vary on a short time scale and without any oeried- 
ici-y, the probaoilistic approach is the only one possible. 
On the other hand, with larger time and space scales, an 
effort in sampling the actual geophysical environment is 
eminently beneficial. In the Sulf Stream region, for 
example, knowledge of the positions of cold and warm eddies 
has proven instrumental for acoustir forecasts. In the case 
of shorter scales, but where oeriodic motions are related to 
a known generation mechanism (e.g., tides) , the same kind of 
hind casting may be possible. 
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It is a*tem pt ed here to shew how 
observed it the Bay cf Biscay may 
ranges, and to assess a possibility 
dence of these ranges or.' relative p 
receiver vs. the internal propagating 
on the direction of. propagation rel 
wave orientation. The internal varia 
assumption is an idealized extract 



gathered from the study of the SAB i 
Therefore, it applies only to this 
ation, and with a wide range of unrer 

This approach may prove feasible 
internal waves can be occasionally 
space using a SA3, or by adequate i 
aues based on a thorough knowledge 
climatology . 

The influence of an internal w 
propagation is especially critical at 
cies, where diffraction and leakage 



the ocean variability 
affect sonar detection 
of analysing the depen- 
ositions of source and 
wave field, as well as 
ative to the internal 
bility taken as a basic 
la t ion of the results 
mage on August 20 1978. 
given particular si tu- 
fa inf y. 

in similar cases where 
sampled, either from 
n situ sampling tech ni- 
of the infernal wave 

ave field on acoustic 
high acoustic freguen- 
are minimal and where 



oure geometric ray tracing can be apolied. The ir.se n if ied 
volume will then drastically vary in che three dimensions. 
Long internal waves will act as long scale time variations 
ir. an observed or predicted SS? (sound speed profile). The 
effect of shoot waves can be approached stochastically and 
lead to establishing statistical conclusions on the range 
variability. The variation in propagation will be dramati- 
cally dependent on ( 1 ) the relative pcsi~ion of source (or 



receiver) with respect to the wave 
direction of sound propagation and w 
and (3) the mean 33? and its variabili 
Two cases are to be considered . 
internal wave field for low freguen 
evaluate variation on long ranges (o 
propagation: in this case the Paraoo 



field, (2) the relative 
ave field orientation, 

Z 7 1. n s C cl. C 3 

(1)A long (e.g. tidal) 
cies (long scales) to 
f order 23Kh) of sound 
lie Equation (PE) hotel 
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has been used. (2) A short 
continental slope) free irt 
cies (short scales) to e 
frequency acoustic wave. I 
run. 

X. IDEALIZED SAVE FIELDS 

The acoustic model wa 
vertical plane the sound s 
profile in a periodic way ( 
was drawn from the XBT da 
smoothed to present a perie 
thermociine and a clinac. 
constant SSP two convergent 
the Mediterranean outflow a 
pressure increase with dep 
made to vary 1 r. depth witho 
ness or its slope (Fig. 3.1 
tier, was modeled as a peri o 
and 2 OKM . 

Increasing this interna 
equivalent to modeling the 
with respect to the wavecre 

sin 
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possib le 
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dimensional 


in Bonification 


1 3 . 179 d 


irections 


from 


diagrams at 


one given 


depth 



(as seen on the 3 AH iraace of the 



err. al wave field a 



— ■ L. ■» 



an treat en- 



valuate interaction with a 



n ich 



n this case the 3EDJSA model was 



s evaluated assuming that in a 
peed profile varied from a mean 
Fig. 3.1) . This assumed profile 
ta, and was then idealized and 
ct surface mixed layer, a sharp 
olcgical deeo profile. For a 
e zones are apparent: one iu = to 

t 1100M, and another due to the 
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- .4 • 


The 


sharp thermociine was 


u* 


any d e 


formation in its thick- 


) i 


the si a 


pe of this depth varia- 


o 

i 


curv e 


with wavelengths of 10 



1 wavelength by a factor "k" is 
ray tracing in a plane oriented 
st s , at an ancle "A” such that 

(A) = 1 /k . 

Iculate transects of the three* 
(or transmission loss) in various 
the source and draw horizontal 
(Fig . 3.3 and 3. h) . 
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TRANS8ISSI 0 N LOSS DSRI7SD FROM THE PR MODEL 



The parabolic equation (PE) model developed by NORDA 
[Brock 1978], modified by NUSC [Lee and Bctseas 1982], and 
adapted (with acre modifications! for use at the Naval 
Postgraduate School, gives a normal mode solution no the 
wave problem in the form of an elliptic equation. The calcu- 
lations made with this model must be considered with some 



caution as 






ste d will 



. — ^ £ 



field results. Principal significance should, therefore, be 
attached to PS modal calculations at the convergence rones 
(CZ) ranges. I 



h — C 3. S 9 


where one source 


was placed over 


rest, no 


3 ^ J ^ ^ ^ £ i. z 


actional deoen— 


and the 


transmission loss 


was equivalent 


const ant 


SS? environment. 


Conversely, in 



the case where the source was placed above an internal wave 
trough, the convergence cones were found to be enhanced by a 
varying SSP. Horizontal diagrams of transmission loss were 
drawn for two levels for the receiver (15 and u 5 ) by inter- 
polation between five vertical transmission loss curves at 
various directions relative to the internal wave field (Fig. 
8 . ' 4 ) . 



C. RAY THACIH G WITH MED OS A 

The MEDUSA Range Dependent Ray Model was developed by 
the Applied Research Laboratories of the University cf Texas 
at Austin under NORDA sponsorship during FY90. MEDUSA uses 
smooth functions (splines) to interpolate between the 
discrete BSP's (given as input) and solves the various ray 
equations numerically for a limited number of character is t ic 
rays (eigenrays) . Eigenra ys are located by interpolation in 
range (not in depth as in other ray tracing models) between 
rays bracketing the receiver in range. The source is 
assumed to lie on the vertical axis of a cviindrical cocrdi- 
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r.at e system, around which the geophysical variables are 
symmetrical. Sigenrays are locatai by ir.tarpolati.ig between 
rays bracketing the receiver in range (not in depth) . 
(Hence any direct modeling of horizontal refraction is 
excluded). The PROF TR? program permits diagnosis of profile 
artifacts produced in the geophysical interpolation routine. 
RATPLT produces ray path plots. 

This ray tracing program was ased in the same way as 
the PS model to compute acoustic transmission in vertical 
planes at various angles with an internal wave field orien- 
tation. Bu- 1 : here, the features which were investigated ver = 
at shorter ranges; the direction dependence exhibited in the 
horizontal diagrams (Fig. 8.3) relates to an eventual 

refraction under the mixed layer, thus cancelling (for rays 
travelling at an angle with the internal wave crests) the 
’’shadow zone". Results are only purely qualitative ar.d would 
need to be investigated from a transmission loss point of 
view . 

Hedusa was used also to input the horizontal sound spee 
profile as observed in the same iiealized infernal wave 
field along an horizontal in the direction of wave propaga- 
tion and at the near, depth cf the thermoclina (Fig. 3.2) . 
The result of ray tracing in the norizontal plane will be 
valid only given the assumption mat the rays do not make 
any depth excursion from the choosen depth of the horizontal 
SSP. This assumption can be made only because the vertical 
ray tracing shows that there is a surface sound speed 
channel, sc that a significant part of the acoustic energy 
will propagate within the range of variation of the mixed 
layer. Above the the fhermccline the sound speed will not 
exhibit any horizontal gradient and no refraction will 
occur. 
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D. GENERAL RESULTS 



Transmission losses (TL) in an idealized internal wavs 
field were found to be equivalent to the TL relative no a 
constant SSP environment if the source was placed over a 
crest. In the case where the source is above a -.rough, and 
for internal wavelengths of 10KM, the transmission loss 
decreases at the convergence zone by IGdS for an insonifica- 
t ion perpendicular to the crests (Pig. 3.4) This anisotropy 
at long ranges, even if not representative of any real ocean 
environment, is at least an illustration that passive detec- 
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Ray tracing shows that there is a possibility of insoni- 
ficatior of the "shadow zone" beneath the theraocline with a 
variable S3? structure. This would affect both passive and 
active search (or escape) tactical patterns a*: lew ranges 

(of the order of 4000*) . Ray tracing in the horizontal 
suggests the possibility of an horizontal duct (particularly 
in deep and narrow internal troughs) , but the assumptions 
made hers are too drastic to allow any conclusion from the 
elementary run made with hedusa. Further work in this line 
should, nonetheless, help to tailor use of any elementary 
(single SS?) ray trace model to interpret some consequences 
of a variable 33? on transmission anisotroov. 
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IX. SAR AS AN INTERNAL OCEAN SAMPLING SYSTEM 



The SEASAT experiment was the first specific cceao.o- 
graphic use of a SAR.frora space. Assessments made a p rior i 
on its monitoring capabilities ha 73 been justified by the 
large number of oceanic internal processes ir observed. In 
less than 42 hours, ir covered an area of 100 million KM 2 of 
the Earth’s surface (35 of them over the oceans), [ Pravdo 



1983 ]. 



The main specific advantages of a SAR as an ocean remote 
sensor system are its all-weather capability (in particular 
with respect to cloud cover, night conditions and sea- 
state), its high resolution and its sensitivity to the 
hydro dynamical surface manifestations of internal activity. 

The cbvicus difficulty encountered with an operational 
use of such a sampling system is the computation time 
required to process the enormous amount of data (some 10 6 



pixels per second). In the case of an orbiting station (as 
proven by -he 3IR-A experiment on the 1982 Shuttle 
missions), the data may be acquired and exploited on board 
and/or kept in memory for a further cround- processing. But, 
for an unmanned satellite, data recording and exploits- ion 
require that an earth station be at radio receiving distance 
from it. 



A. G SHE HAL RE QUIRE M ENTS PO R A SAR DCEAN MONITORING SYSTEM 
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It just y i 7 e a con tin uous radar coverage over 1 sr gth 
scales greater ‘•.han the mesoscale features than must be 
detected; This a 3AR designed solely for sampling the swell 
may well work by sampling 10KM squares of ocean each 100KH 
but will be unable to observe internal wave fields. 

Another quite obvious characteristic is related to the 
process cf creating a synthetic aperture. The orbit charac- 
teristic of the spaceborne vector must provide a relative 
velocity with respect co the earth's surface, thus a 
geosynchronous satellite cannot be used. Nevertheless 
repeating orbits or qua si-geosynchronous orbits ( i. e. , 
near-polar orbits with periods equal to the earth's rota- 
tion) could combine the convenience of a specific coverage 
over a given meridian area and the necessity of a relative 
sd ee d component to permit use of the Doppler integration 
scheme [Tomiyasu 1982]. ? or the investigation of internal 

features, a sun-synchronous orbit would no 4 - be optimal 
because a given area would always be monitored at the same 
time each day. 

The radar frequencies at which the surface wave field 
responsible for the 3ragg effect will effectively interac” 
with internal features may vary from the S-band (a few 
centimeters wavelengths) to the L-band (a few decimeters) . 
High frequencies should yield better wind-speed sensitivity, 
but would increase the difficulty of inf er oretat ion since 
the dynamics of gravity-capillary waves are more complicated 
than those of pure gravity waves. 

The viewing angle may itself vary from at least 15° to 
60° (SIR-3, Spaceborne Imaging Radar of fh- coming Space 

Shuttle flights) , will have a capability to vary mechani- 
cally the viewing angle within nose limits) . If seems 
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13 SITU CO MPL2HENTAR Y 3 INOLTANEOO S PROBING 



To analyse c haract eri sties cf she internal features 
displayed by a 5AR, the mean density structure mus~ 'os 
known. It has been shown (Cf chapter VI) that the nixed 
layer depth, the amplitudes of high frequency internal 
waves, and the location of fronts can be inferred from 
consequent surface patterns. If would therefore be possible 
to conduct such an analysis in real time if some X3T 
sampling was made nearly syncpticaily and could be blended 
with climatological data. 

Obviously time series of temperature structure would be 
more efficient still by giving synoptic information on the 
frequency spectrum and vertical modes of the internal 
motions. A renote, drifting thermistor array as described in 
[N03C 1930], if launched in an appropriate position with 
respect tc the a pri ori internal wave field, could give both 
ar. absolute reference to radar renote sensing results and 
fill in the gaps between two consecutive passages of the 
satellite. (A radar trar.spondor could add, at lew costs, a 
geographic reference.) 



C. OPERATIONAL PROCESSING AND EXPLOITATION 



The analysis of a SAP image may be made 
levels corresponding to different goals. 



on two differ enc 
time constraints 



and processing facilities (Fia. 9.li . 
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cion by meteorologists of a visible or 
image cf the atmosphere, would probabl 
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similar features is capaple of giving 
relative properties cf a given internal 
the geographic misplacement of fronts, 
measurements of interpack.ee separaoion. 
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variation of wavelength, couli be is ad ir. simple araiy 
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probably help to tailor processing algorithms, eventually 
using the raw data directly. Once the features of interest 
are defined by their approximate locations, dimensions and 
average characteristics, the regained data could be 
extracted as a subset and processed according to specific 
methods. For example, if a front is known to move parallel 
to itself and to be charar terizsd on a SAP. image by a band 
of high radar intensity due to convergent enhancment of the 
surface waves, the mean radar intensity averaged parallel to 
the front would be encuoh to locate its oosition (and inten- 
sity) . If the presence of high frequency internal waves is 
assumed, a procedure to extract their mean direction and the 
distribution of their wavelengths would be to analyse a set 



of pixels over a limited area by 
Fourier Transform technique. 



a two-dimensional 



Fast 
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X. SaHMARY AND RECOMMENDATIONS 



The SAR image gathered by SEA3AT over the 3ay of 3 i sc ay on 
August 20 1978 showed a number of internal features which 

were generally in good agreement with analysed synoptic and 
climatological data and could have alternatively be a n used 
to infer valuable sea -truth JLn operational, real-time 



sampl ing. 

Wave packets propagating offshore from the shelf break 
had an inter- packet spacing on tie order of 23KY corre- 
sponding to a semi-diurnal generation at the main canyons 
indenting the shelfbreak at low oarotropic tide (maximum 
ebbing surface flow) . Speeds computed from observed density 
profiles were 0.55.1/S +/- 0.5. over the slope, which 

compares well with the wave speed of 0.511/5 calculated for 
those packets from the SAR image. These high frequency 
internal waves have a period of the order of 701TN. 
Assuming those waves became linear far enough from their 
generation points, a mean depth of the mixed layer of 201 
was inferred and spatial variations ascertained from clima- 
tological data and observed group spacing. 

Drastic spatial variations in tie wavelength observed in 
a u i ven gro up and between differs nt ones were a r a 1 v sed , and 
assumed to be due to current shear and horizontal density 
gradients. They coincided with a sharp thermal front at the 
shelfbreak and an alongshore surface geostrophic jet. A 
line of surface enhancement of 3ragu waves corresocr.dsd to 
this front, which was analysed by quasi-synoptic X3T 
sampling. This front is generally present in summer condi- 
tions when a rorresocndinc UDwelliac breaks the surface. 



23 



Nonlinear descriptions of the internal waves in each 
internal wave packet yielded estimates or their amplitudes, 
probably as large as 15M on the continental slope. 

A tentative application to underwater acoustics shows. 



under severe simplifying as sum 
in such a situation -are as sen 
that sonar search patterns (or 
routes) could be greatly enhar. 
could be sampLed on a regular 
Adequate operational X3T samp 
information on this variabilit 
areas where internal waves are 
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corresponding to half a period 
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expected. Internal wave cli 
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ptions, taat detection ranges 
daily direction-dependent and 
, conversely, submarine escape 
ced if these internal features 
basis and taken into account, 
ling could give seme in -situ 
y. It is recommended that, in 
known to occur, consecutive 
time (or distance) interval, 
(wavelength) , thus giving an 
1 wave amplitude and ccr.se- 
deteotion ranges that can be 
matology developed using SAB 
igning a sampling plan. 



The following ocean features of operational interest can 
be inferred from the SAR image studied here: 

a) An area of intense internal activity is clearly visible 
on the continental slope in the form of short internal wave 
groups propagating from the heads of the major canyons and 
generated approximately at the time of lew fide at Brest. 
These perturbations probably correspond with a Larger scale 
widening (weakening) of the thermocline and periodic mixing 
processes. The possible operational consequences are 
(1) drastic changes ir. detection ranges due to a varying 
thermal structure in time and space and a probable increase 
in ambient noise level due to bottom stress and turbidity 
currents, and (2) sudden mo dif ications in buoyancy that may 
affect underwater vehicles (e.g., submarines, torpedoes) . 

b) A well-defined front, located at the upper mart of the 
shelf break, is ao oarer.. on the BAR inace as a stror.c 
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localised enhancement of the high frequency waves. It corre- 
sponds to the well-known shelf break upweliing in the summer 
in the area, which has been demonstrated by synoptic X3T 
sampling. Insonif icatior. through such a front is essen- 
tially dependant on the relative positions of source and 
receiver. 

c) In the open sea, several fields of internal waves are 
organised in groups spreading from the shelf. The general 
direction of the crests is parallel to the isobaths and 
creates anisotropies in sound speed structure and insonifi- 
cation volume. 

d) Other features, not considered here, have eventual 

tactical importance: (1)the localisation of areas of calms 

and therefore mirror sea surface for visual search in ASH 
and SAB (Search and Rescue) missions and site selection for 
weapon launch, and (2) the coastal refraction of the swell 
giving a possible surf hindcasting capability for mine- 
sweeping and amphibious landing operations. 

The coming SIB-3 missions will cover the Bay of Biscay 
and consequently provide some information on the time vari- 
ability or the characterist ics appearing in this single SAR 
image. Combined international effort to increase synoptic 
in- situ data acquisition at this time should be eminently 
ben ef icial. 
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dramatically change the parameters of undersea warfare. 
Tver y effort should therefore be made in allied navies to 
understand and master these consequences, from both scien- 
tific and operational points of view, to make sure that no 
dangerous and destabili zing consequences could result from 
their universal exploitation. 
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